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Page 10, Section 2.1.1 First paragraph: change date to data. 
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Scattering and Diffraction 

Page 10, Section 2.1.1 First paragra^ : change date to data. 

Page 59, Equation 14. Ad^.: for n<m 

Page 60, Equation 16, first part: Add for n<m. 

Equation 16, change to: 

Page 63, Fourth line: Change to read as follows: 

as a ficticious unit area surface normal to the incoming or out- 
page 67 , Paragraph 6.4.1 line 8, change to read as follows: 

too cumbersome in the analysis cf the fine-scale diffraction 
from the struts 

Page 90, first line. Change 6=0, to read 6^^0 
Page 92, Equation 74. Change to: 

(p = cos ^ { tanwtan (1^/2) } (74) 

Page 100, 5% Diffuse at 88°: change from 3.5E-11 to 1.5E-11 

Corresponding percentages now are: % Diffraction 99. 

% Scatter 1. 
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Page 102, Second line. Change % Diffuse to 5-o Diffuse 

5"u diffuse at 24® » clmnge from 9.5E-4 to 9.5Ep5. 

Corresponding Percentages now are % diffraction 1. 

% scatter 99. 

5% diffuse at 30®, change from 9.8E-E/ to 9.8E-6 
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Corresponding percentages now are; % diffraction 4. 

% scatter 96. 

Second line from bottom, change % Diffuse to % Diffraction 
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Corresponding percentages now are: % d ffraction 9. 

% scatter 91. 

Page 102, Second summation with Martin Black at 88® 
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Page B-10 Change line 15 to: 

(x.y) = (0 ± R). 

Page B-11 Change asymptotcc to asymptotic. 

Page B-13 First equation (the one before B-46) . The third exponential 

term in the brackets should have a minus (-) in the exponent 
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1.0 INTRODUCTION 

This report contains the results of stray-light transmittance analysis 
on Perkin-Elmer* s design of the Infrared Astronomical Satellite (IRAS), 
design mir;;bcr 693"10000, Revision E, dated May 9, 1978. The system was 
evaluated for scattered radiation propagation with the use of the APART 
(Arizona's Paraxial Analysis of Radiation Transfer, version 6), and the 
propagation of diffracted energy with both PAOE (Paraxial Analysis of 
Diffracted Energy) and Perkin-Elmer' s GUERAP II programs. The results of 
scatter and diffraction are first presented separately, with the combined 
transmittance values being presented later. 

The scattered radiation aiialysls was performed by using both a 5% 
diffuse black on all the baffle and vane surfaces and by using a math- 

i- 

ematical model of the Martin Black scattering characteristics. The results 
of the IR scattered- rad iat ion analysis show that the majority of radiation 
comes from the inner-secondary baffle and the object side of the aperture 
stop. In the visible wavelengths, the primary and secondary mirrors are 
the dominant sources of unwonted radiation at all off-axis angles. 

For all wavelengths, diffraction effects are dominant only at the 
large off-axis angles, except for a few specific casas at smaller angles. 

This is somewhat contrary to popular opinion which holds that at the 
long wavelengths diffraction effects will predominate. This report 
clearly shows that at the longer wavelengths the diffraction contributions 
go up significantly; however, surface-scattering characteristics are also 
larger, resulting in more unwanted energy reaching the Image plane due to 
scattered radiation. It is the comparative increase which determines 
which propagation process predominates. 


I 
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The followinq tnsks were outlined In the Statement of Work;' 

1. Perform an lnd<*oendent analysis of the off-axis rejection of 
the IRAS Telescope System (including sunshade, optical sub- 
system, and field optics) considering the effects of both 
scattering and diffraction over the operating wavelength range 
of the telescope. The telescope design analyzed shall represent 
the flight design to the maximum extent possible. The optical 
subsystem^ Is defined to Include the telescope optics, struc- 
ture, and baffles. The requirements for off-axis rejection 
shall be as defined in Specification 2-26iil2 "Performance 
Requirements for an Infrared Telescope System for the Infrared 
Astronomical Satellite (IRAS)" Revision 5, dated September 15. 
1977 or the latest modification thereto. The analysis shall 
include but not be limited to the following elements: 

(a) Computation of the off-axis rejection of the optical 
subsystem by itself for direct comparison with the 
Perkin-Elmer ariclysis. 

(b) Computation of the azimuthal variation of diffraction of 
the telescope for comparison as in (a) above. 

(c) Computation of the effects of the field optics Including 
both the cavity behind the field lens and the aperture 
that proceeds the detector cavity. 

(d) Computation of the system off-axis performance at 
632.8 nm and 10. nm. 
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The analysis In each task shows, in most cases, c]Ood ocjreoPH'nt with 
Perkin-Elmer ' s original analysis, even though the systems analyzed were 
s 1 ight ly di f ferent. 

The specular sunshleld effectively blocks off the solar rndiaiion 
for off-axis angles greater than However, because of its specular 

coating. It may collect and reflect unwanted radiation from near off- 
axis sources into the system. The field masks and lenses in ttu' final 
flight system will alter the propagation piths to the individual dotectof 
In all cases this will result in better off-axis rejection. For a few 
off-axis positions, the nvask;. will provide the crucial reductioti required 
to meet the specifications. 

The stray light requirements used in this analysis are those defined 
In Spec i f Icat Ion 2-261412, which are repeated here; 


3 • 2 . 2 . 3 Stray Light Rejection 

The Telescope System stray light rejection require 
ments arc defined as follows: 


a. Let P{0) be the power (watts) from an 
unwanted point source detected when the 
Telescope System's line of sight is 
displaced an angle 0 from the point source. 
P(o * 0) * P(0) is then the power that 
could be detected if the point source 
were imaged directly on the detector. 

b. The normalized off-axis attenuation A(o) 

is defined as ** 


A(0) 


Q 


P(0) 


sr 


1 


where is the solid angle (sr) subtended 
by the detector. The Telescope System shall 
have a normalized off-axis attenuation equal 
to or less than the values of A(G) tabulated 
be 1 ow in Table 3. 2.2. 3. 


k 


Table $.2.2.3 
REQUIRED A(0) sr‘* 

Spectral Band (pm) 


G 

0.1* - 0.9 

P - 15 

0 

1 

kb - 81 

87 - 

118 


1 X 10 

7x10° 

3 X 10-1 

k X 10-2 

2 X 

10“^ 

2/«“ 

3 X 10-1 

1 X 10-5 

8 X 10-° 

X 10-° 

5 X 

10-° 

60“ 

3 X 10-7 

5 X 10-° 

1 X 10-7 

2 X 10-7 

1 X 

10-7 

88“ 

1 X 10-7 

2 X 10-° 

X 

o 

1 

U X IO-° 

4 X 

10-° 


In any computerized stray-light analysis, the model of the system Is 
very Important, Including the scattering characteristics of the surfaces. 
The surface elements that were used in this analysis are shown In Figs. 

I to 3. Figure I shows those elements that were tnodeled as existing In 
object space (space one). Table 1 describes each surface according to 
the numbers that appear In the three figures. Because of the nature of 
the APART program, the appropriate elements must be entered Into each 
space of the system. As can be seen In the three figures, the elements 
In each space nr i not be the san'e. 

Elements numbered 5 and 6 are the sections of the primary baffle 
which have vanes. In APART, the locus of vane tips Is entered as the 
"surface" of the element. This accounts for the odd shape of sections 5 
and 6 relative to section 8. 

Section 8 has some localized vane structure to shield rivets along 
the barrel baffle seam that Is not accounted for In the analysis. These 
structures are comparatively small and should not adversely affect the 
resul ts. 





12jt< 
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Table 1. Description of Objects Used In APART, 


Object ^ 

1 

2 
3 

5 

6 

7 

8 
9 

10 

n 

12 

13 

lA 

15 

16 

17 

18 

19 

20 


22 

23 

2h 

25 

26 

27 

28 


« o :rm at ■ar,m «r.» <i «'Ui^ar -V » s. «, • k <«' «C'1K3«.9c.T’< 


Poser Ipt Ion 

Elliptical tilted diffracting edge 

Sliced at 31" conlr*;l reflecting shield 

Disk as rcullotors 

Entrance port diffracting edge 

Front section of njain tube v^/here first 5 
vanes are 88.5 separation 

Second front section of main tube where 8 
vanes are spaced AA.2 tnm apart 

Right side of the last vane of object 6 

Vancless side wall of the meiln baffle 

Vane at right side of object 8 

Small right side of object 9 

Baffle extending from aperture stop backward 

Left side of aperture stop vane 

Diffracting aperture stop tip 

Right side of aperture stop 

Small cylinder extending from aperture 
stop to primary mirror 

Cylindrical outside of secondary structure 

Conical outside of the secondary baffle 

Secondary baffle tip 

Inside of secondary baffle 

Small extension from secondary baffle 
toward secondary mirror 

Back side of spider support 

Outside of Inner conical baffle 

Inner obscuration vane at the primary mirror 

Diffracting edge at tip of Inner conical baffle 

Diffracting edge at tip of inner conical baffle 

Cylindrical end of inner conical 

Middle part of the inner conical baffle 

Cylindrical right end of inner conical baffle 


Sect Ions 
I X 5 
5 

<t 

1 

5 

k 

2 

5 


10 x 3 
5 X 5 


it 

2 
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Table I. 


Description of Objects Used In APART, Cont. 


.TiS- ■» » dK -a 




Object # Poser I pt ion Sections 


29 

Steep cor»e at right end of Inner conical baffle 

2 

30 

Image baffle 

2 

31 

Image 

5 

32 

Secondary mirror 

<i 

33 

Primary mirror 

6 

3'* 

Not used 


35 

Primary mirror as used In space three 

6 


6 



2.0 ANALYSIS OF THE OPTICAL SUBSYSTEM (OSS) 

2.1 COMPARISON WITH PERKIN-ELHER' S ANALYSIS 
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2.1.1. The Visible Band . The analysis of the optical subsystem 
(OSS) was performed to permit a comparison with the Perkin-Elmer analysis 
as reported In PE Report No. 13616. Initial evaluation of this report 
and a listing of the input data decks revealed some conflicting informa- 
tion. Figures I and 2 from the report are reproduced here as Figs. 
and 5. By observat-lon, one can sec that the size of the Inner-primary 
baffle near the primary Is much too large and the aperture stop appears 
to be too small. A check of the data input listing, Fig. 6, shows that 
neither is the case. As was explained In the Introduction, object 2 In 
the PE report should have the locus of the vane tips used as the "surface”. 
The drawing is inconsistent with the actual (and correct) date Input. The 
only error in PE's input data is the reference values of the marginal and 
chief rays for each of the spaces. The detailed effect of this error on 
the calculations Is quite complicated; briefly. It will shift the loca- 
tion of the image of the objects. 

Table 2 shows the percent of power contributed to the full detector 
area (r ■ 5.04 cm) from each of the individual objects. The black coat- 
ings were assumed to be 5? Lambertian diffuse. Ths mirror coating had 
a BRDF of I.75E-I sr"^ at 6 - 6^ * .01 with a 6”^ falloff. This value 
Is loo times higher than the BRDF specified by PE for lO p. However, 

PE later^ did specify a (X/X )“*^ scaling at X **10 p. The mirror BRDF 

o o 

used is thus about 2.3 times lower than the appropriately scaled data. 
Subsequent analysis reported here will Include the BRDF scaled according 



(Units Meters), PE Data. 




. ^E-RKIN ELMER 

Report Mo, 



Fig. 5. Model Used for APART Program Calculations, PE Data. 
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6. Input Data Listing from Perkln-Elmer's Analysis. 


Table 2 . Z Power Contributions to Image for Visible Wavelengths and 5 Z Diffuse. 
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to the PE formula. In either cat»e, the significant contributors are the 
Inner- secondary baffle, the prinwry, and the aperture stop. 

A detailed critique of the design or an analysis for each off-axis 
position will be avoided here as this analysis Is to compare results with 
those obtained by PE. Significant differences In the principal propaga- 
tion paths will be highlighted. When the .'light version (OSS plus sun- 
shield, field optics, and masks) is evaluated, there will be greater 
detail so that the effects of suggested design changes will be readily 
recognized. 

The above analysis, with a 5"^ diffuse black coating, was repeated 
using a niodjl of Martin Black as the surface scattcrer. A brief descrip- 
tion of the nxjdcl is given In Appendix A. This model accounts for the 
higher forward scatter and lower back sea; ler that Is characteristic of 
Martin Black. The results with the Martin Black model (Table 3) show 
some significant changes In the propagation paths duo to the above char- 
acteristics. The A(0) values are plotted In Fig. 7 f-''r both runs along 
with the specification for the 0.1* to 0.9 u ban.l. 

A comparison of PE results (Fig. 8) with our analysis shows excellent 
agreemeut except for one data point at 10®. PE reports that the major 
contributor at this angle should be the inside of the secondary baffle. 
However, this was added as a separate hand calculation.^ APART calculated 
the power focused onto the Inner-secondary baffle; It then computed the 
radiation scattered forward toward the Image of the detector as Imaged 
by the secondary and also the energy backscattered to the primary and 
through the secondary before reaching the Image. The area directly 
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seen from the Image did not receive direct Illumination, thus, the hand 
calculation appears to be in error. 

Figure 9 is a very simplified drawing of the rays that can reach 
the secondary baffle. The main baffle tube blocks out all the radiation 
to the sections on the secondary baffle that can be seen directly. 



Fig. 9. The Left Sections of the Secondary Baffle do not 
Receive Direct Illumination. 

There is one last point to discuss in the comparative analysis of 
the visible band results. This APART analysis predicts the contributions 
from the front side and back side of the aperture stop is the reverse of 
that presented by PE. The magnitude of the power contributed also appears 
to be reversed. Only a small portion of the back surface can receive 
radiation compared to the entire front surface, which is seen directly. 
With Lambertian coatings, then, the front side must contribute more. 

2.1.2. The 8-l5u IR Band. The OSS was evaluated for its off-axis 


rejection using a Lambertian diffuse coating. The mirror coating had 
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.1 BROF of 1.75 E-3 at • .01 with falloff. Is the sine of 
the angle of tncldcncc. while B Is the sine of the angle from the surface 
normal to the observation point. Table ^ shows the contributions from 
each object. One can readily see that the mirror scatter plays a second* 
ary rold In the IR. In this analysis. It Is by definition (and generally 
accepted as a true statement) that mirror scatter goes down with Increas* 
Ing wavelength. Therefore, unless the BROF of the black surfaces also 
decrease with wavelength, one should not expect significantly better off- 
axis rejection at the longer wavelengths. 

The major contributors are the Inner- secondary baffle, the aperture 
stop, and the Inner-primary conical baffle. Note that the back side of 
the strut contributes only 5*6^ of the total energy and that Is only at 
one angle. The or Icntat Ion of the struts Is shown In Figs. 10 and II. 

In particular, one of the struts Is aligned with the peak of the sun- 
shield tip. The values calculated In this analysis have the strut and 
the off-axis point source In the same plane. For small off-axis angles 
the unwanted radiation is (nearly) focused onto the strut, making this 
azimuth the worst case. 

Table 5 shows the same analysis but with Martin Black on the baffle 
surfaces. As was the case with the analysis In the visible band, the 
percentage numbers vary but the same objects are critical to the system’s 
performance. This Is due to the variation In the forward and backward 
scatter characteristics. In most cases this shifts the percent contri- 
buted from a back scatter path to the forward scatter path. Figure 12 
shows the relative A(0) values for the two runs along with the specifica- 
tion. The Martin Black A(0) values are more often below the IRAS 8-I5u 


spec line. 
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Fig. 10. Strut Orientation of the IRAS System. 
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Model of IRAS Used in the APART Calculations with Vanes and Profile of Sunshade. 
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Fig. 12. A(0) for 8 - 15 p Band of OSS. 
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Figure 13 compares the PE results with the 51 results. The contri- 
bution calculated by PE, as Indicated by the dashed line at 10 degrees. 

Is much higher. The contribution from the front side of the aperture 
Is about 10% higher than predicted by PE for the back side. These are 
the same values as calculated In the visible, because the 5 % black 
hasn't changed with wavelength. 

2.1.3. Black Coatings at Long Wavelengths . When the Martin Black 
coating was used evaluate the OSS In the 8-l5u band, the predicted 
values were. In several cases, above the 8-l5w spec line. At the longer 
wavelengths the spec line moves down to lower rejection values. However, 
existing data indicates that the hemispherical diffuse reflectivity In- 
creases with wavelength (Fig. \^), This implies that at some off-axis 
angles A(g) will be above the spec lines for all the IR bands. 

The above staten»ent must be tempered with the following additional 

o 

statements: Other measured data exists which indicates that the increase 

is not as pronounced (Fig. 15). At the longer wavelengths the surface 

roughness relative to the wavelength, Is much less. As with the mirrors 

$ 

one might suspect that the diffuse BRDF should drop significantly while 
the specular component will increase. It Is not known whether the data 
In Fig. I^» Includes both the specular and diffuse component. 

In any case, it is imperative to have measured BRDF data as a func- 
tion of the input and output angles. Such data is not presently avail- 
able for the wavelengths above 10. 6p. Without it, the validity of the 
evaluations at the longer wavelengths is questionable. 
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Hemispherical Diffuse Measurements as a Function of X. 
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In the analysis of the flight system, discussed later in this report, 
it will be assumed that the BRDF profiles at ail wavelengths remain the 
same, only the magnitude varies with wavelength as indicated in fig. U. 

If, as discussed above, the coating exhibits a marked specular 
reflection at the longer wavelengths, this could significantly alter the 
predicted attenuation factors. Most likely the specular component will 
not be a delta function but exhibit a relatively high BROF over a fairly 
targe solid angle (» .02 sr*M> Under such circumstances, the forward 
scatter from the secondary baffle towards the secondary mirror would bo 
significantly higher and specular reflections from the inner conical 
baffle towards the field mask would cause additional problems for large 
off-axis angles. It would then be desirable to eliminate all possible 
near-specular paths. 

3.0 SUNSHIELD ANALYSIS 

This portion of the analysis was performed in two stages, in this 
section the characteristics of the sunshield itself wilt be explored. 

The overall effect of the sunshield will be discussed later as part of 
the evaluation of the overall flight system design. 

The sunshield design can be seen In Fig. 10. Figure 16 shows a 
profile of the shield. The design of the shield has been detailed else- 
where;*® basically It blocks sunlight for off-axis angles greater than 
60® when the sun Is in the plane of the tip of the cone. Earth light 
is rejected at 88® for off-axis angles In the opposite direction. 

The underside of the sunshield has a specular surface to reduce the 
thermal loading on the telescope. When the unwanted source is more than 
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off axis only diffuse scatter and self-emitted radiation will enter 
the main tube, for specular surfaces* the emissivity is low and the 
diffuse scatter Is low, hence the heat load will be low and the attenua- 
tion high. 

However, at small off-axis angles unwanted light is specularly 
reflected Into the main tube. Increasing the heat load and the A(o) 
values. The problem Is three-dimensional , but Fig. 16 will help to 
clarify the problem. Ray I from a source 5 . 7 ® off-axis will Just reach 
the entrance part of the main tube. So the Increase to the specular 
reflected energy should start at guite small off axis angles. Ray 2 is 
from a source 23. off-axis. Beyond this angle, specular paths do not 
enter the entrance port and the Input energy will fall off. 

The three dimensional nature of the problem also causes skew rays 
to be focused Into the system. To analyze the problem quickly, several 
scale models of the sunshleld were made using specularly coated Mylar. 
Transmittance measurements were made as a function of the off-axis posi- 
tion of the source (Fig. I7). The peak Input power is for a source point 
25® off axis as shown In Table 6 . This Is very unfortunate because the 
A(o) for the OSS Is Just slightly above the specification for off 
axis angles of 2l«-60* when the system was evaluated using Martin Black. 
The sunshleld will cause the values to go even higher. 

By changing the angle of the sunshleld, the peak value due to spec- 
ular reflection may be moved to a less obnoxious off-axis angle. Depen- 
ding on how this Is done, It will usually affect moments of Inertia and/ 
or the radiators that presently exist. 




Measured Perfornance of the SonshteTd HodeJ. 
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Table 6. Data from Measurentcnts on a Sample Sunshade. 


Off-axis 


Ang le 

£ 

0 

1.00 

5 

1.03 

10 

I.IO 

15 

1.29 

20 

1.59 

25 

1.69 

30 

1.64 

35 

1.51 

40 

l.40 

45 

1.18 

50 

.99 


♦in ■ »o ‘ 

di,*, Is the power Into the main baffle sections 
IN 

is the power into the system without the sunshicid 


4.0 FIELD MASK AND LENSES 

To analyze the effects of the field stop and the field optics, a 
single detector was selected and located on axis. The field of view 
from the detector was calculated using ACCOS by running rays from the 
edge of the detector surface backwards through the lens, and using the 
field stop as the limiting aperture for the rays, determining the direc- 
tion they leave the field-optics set. This was done for the band 1 
field-lens assembly with the data as supplied by Ball Brothers Research 
Corporation in a 16 May 1978 letter number 86563.78.0.0038. A picture 
of this process is shown in Fig. 18. From this data, it was determined 
that the transfers from the objects along the inside of the inner primary- 
mirror baffle could not reach the detector (Objects 26, 27, 28 , and 30). 
These transfers were eliminated from the analysis for the <R single detector. 




View from the Detector through 
the Field Mask and Optics. 



The field stop is oversized allowing all the stray radiation from the 
secondary mirror (object 32) and objects imaged through it to reach the 
detector. These transfers were changed to account for the smaller size 
detector rather than the whole Image surface. The limiting field of 
view from the detector casts an el i ipticai-shaped hole onto the inner- 
secondary baffle. To model this transfer, the size of the eiilptical 

I 

aperture projected o.jit to the secondary baffle was determined and entered 
as an obscuration in APART, with the detector being modeled as a roctangu- 
lar disk. 

k.l DIFFRACTION EFFECTS 

The size of the field stops must be oversized because they are in 
the far-fleld diffraction region of the system aperture stop. The field 

i • 

stops are, however, not in the far-fieid region for power transfers coming 
directly from the inner-secondary baffle, where here we must consider near- 
fieid diffraction from a linear edge. In this case, the light rays dif- 
fracted around the edge will be of secondary Importance compared to the 
direct rays near the edge which will pass undiffracted to the detector. 

The inner-secondary baffle Is the only object in which edge diffraction 
would produce any measurable contribution. 

5.0 ANALYSIS OF THE FLIGHT DESIGN 

The preliminaries of the scattered-light analysis on the flight 
design have been considered. What remains is the scatter and diffraction 
analysis on the flight-design system as a whole. This section contains 
the stray-light analysis of the flight system including the OSS, sunshield, 

. f 

field mask and lenses. The diffraction analysis will be presented in the 


next section. 
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5. 1 SYSTEM MODEL 

The physical size and shape of the objects are shown in Figs, i, 2, 
and 3i except for the detector area which is shown in Fig. 18. 

The BRDF values of the mirrors are according to PE's formula as 
described in Appendix A. By definition, there is no wavelength scaling 
with the S% diffuse surface. When Martin Black is used the values and 
profiles remain constant except for the 100 micron band. In this band 
the DRDF profile is ASSUMED to remain the same while the magnitude is 
increased by a factor of 10. This results in an inc«*ease In the A(0) of 
10^ where N is the number of Martin Black-surface scatters encountered. 

The scattered light from the sunshleld becomes a factor only for 
very large off-axis angles (s88®). The specular component Is of no 
consequence at these angles. No measured BRDF values were available for 
this surface, therefore it was assigned a Lambertian diffuse component of 
.001 for all wavelengths. At these angles the sunshleld is the only sur- 
face that is directly illuminated so the A{@) values will scale directly 
with the sunshlelds surface scatter. The .001 BRDF yields A(0) values 
well below the spec line and, as will be seen, diffraction effects pre- 
dominate at these large off-axis angles. 

5.2 SCATTERED LIGHT ANALYSIS Q.k-O.S BAND 

The contribution of power from each object to the image plane Is 
shown in Tables 7 and 8, for the 5% and Martin Black coatings respec- 
tively. 

At 5® the direct scatter from the directly illuminated primary and 
secondary mirrors are the major sources of scattered light. 
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Table 8 . t Table for 0 . 4 - 0 . 9 u Band with Martin Black for' Flight Systaa. 
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At 10*, th« Mcondary it no longer directly illufliineted. The prlmery 
mirror end Inner-tecondery baffle ere now the major contributors of 
scattered light. The mirror scatter remains fixed for both the S% diffuse 
and Martin Black calculation, so the relative Importance shown In the 
percent table Is due to the change In the black-scattering characteristics. 

There are three paths of scatter from the secondary to the detectors 
(Fig. 19); 

1. Directly from the left most sections (20t) of the baffle 
to the detector. 

2. To the image of the detector as reflected in secondary (right 
most sect Ions - 20%) 

3. To the Image of the detector as reflected by the primary 
and secondary (right most sections - 20%). 

The projected solid angle of the detector, as seen from the source 
sections, is 100 times less from the last two positions when compared ’ 
to the directly seen area. However, the left most sections are not di- 
rectly Illuminated at, or beyond, 10 degrees. 

• ■* 

With the 5% black the forward and backscatter contributions are 
approximately equal. This is unrealistic. With the Martin Black mode-l 
the black-scatter path (Path 3) drops by a factor of about '0 while the 
forward scatter BRDF goes up by about 100. The forward -scat ter path is 
the dominant propagation path because of the near specular forward scat- 
tering c laracter istics. 

For angles between 17® and the primary mirror and aperture stop 
are both Illuminated. The primary receives considerably more power and 
is thus the major contributor; its relative contribution decreasing as 
the entrance aperture shades it at larger angles. 



For «ngl«t larger than 2^* no objaett art i I lumlnatad that art 
taan from tha Imaga. Up to 88* tha main tuba (objaett 5, 6. 8) ara 
Mlumtnatad and thata In turn scatter to all of tha objects seen by tha 
visible detector. At tha 30* angle object 8 can transfer power directly 
to the inner -secondary baffle which causes Its contribution to be higher 
than for larger angles. For all angles larger than 30; the contributions 
remain fairly constant with only the magnitude dropping with higher angle. 

At angles greater than 88* tha specular heat shield is illuminated 
(assumed reflectivity p ■ 0.001) and this drops the scatter another k 
orders of magnitude. 

In summary, the inner-secondary baffle, inner conical baffle, and the 

mirrors are the dominant sources of scatter. The BROF at 6-B^ " 0.01 of 

0 

the mirror is Q.k]k sr”* which Is high even for average quality mirrors. 
Somewhat better performance can be expected. 

Figure 20a shows the plot of A(0) for the 5% diffuse, Martin Black 
model, and the spec line. The only problem in meeting the specification is 
at the 60 * angle; diffraction is always below the scatter. 

5.3 SCATTERED LIGHT ANALYSIS (ilu BAND) 

For ail the IR bands, the importance of mirror scatter drops signifi- 
cantly. In Ideally baffled systems, the A(0) values should be determined 
solely by the optical surfaces (mirror scatter) or by diffraction. However, 
in this system the limiting factor Is either the black surfaces or diffrac- 
tion. Based upon the available information, the hemispherical -diffuse 
scatter for Martin Black remains relatively constant until the lOOu IR band. 
The A(0) value will not drop significantly If the principle contributor is 
a black surface. 


kk 
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Fig. 20a. A(9) for Martin Black and S% Coatings in the 0.i(*0.7y Band. 
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S.3«l. 5* •nd 10* Off-axU Source Potitlont . Th« principle peih 

of propegeClon Is from the source to the inner-secondery baffle (after 
reflecting off the primary) to the detector as seen In reflection. This 
Is the same path which was discussed In detail In the discussion of the 
visible band. 

5.3.2. 17*. 2k* and 30* Off-axis Source Positions . The dominant 
path of scatter Is from the front side of the aperture to the detector. 
This Is a bacRscatter direction for near-normal angle of Incidence; a 
near optimum condition for Hart In Black. The path Is the result of the 
stop location. Because the stop Is not at the secondary, the physical 
size of the secondary mirror must be oversized to accomodate a field of 
view, allowing out-of-field elements to be seen from the detector. 

At 21)” and 30* off-axis angles, power is directly loaded onto the 
rear sections of the main tube (elements 6 and 8). From there energy 
scatters directly to the inner- secondary baffle. Half of the radiation on 
the secondary baffle comes from the last section of vanes on object 6, 
while the other half comes from object 8. The radiation on the inner- 
secondary baffle then scatters towards the secondary mirror and reflects 

to the detectors. This second-order path is significant enough to almost 
equal the first-order path from the aperture stop to the detectors. 

The specular reflections off the sunshieid have their greatest 
Impact at these angles and continuing up to about 50* off-axis. The 
amount of increase is relative to the values shown in Fig. 19. 

5.3.3. 60* to lOc Off-axis Source Positions . At 60* off-axis the 

propagation path is from the source to the front section of the outer- 
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primary baffle, to the Inner-secondary baffle (via a reflection off the 
primary) and then forward scattering to the detector. The specular reflec- 
tions of the sunshield are no longer significant because they do not enter 
the main tube. 

For off-axis angles greater than 88”, the first collecting element 
must be the specular sunshield. This radiation is scattered ic the vanes 
on the main baffle, then to the Inner- secondary baffle, and finally to 
the detector. • 

S.4 SUMMARY OF THE SCATTERED- LIGHT ANALYSES AND RECOMMENDATIONS. 

Figures 20b to 20e show the predicted performance of the IRAS system 
in each of the wavebands. From the previous discussion and the percent 
table (Table 9) for the 11 micron band, one sees that the mirror scatter 
Is not significant. To reduce the scattered light, one or more of the 
following steps must be made: 

1. Find a better black coating. 

2. Change the projected solid angle between the inner-secondary 
baffle and detectors and also between the aperture stop and 
the detectors. 

3. Reduce the power that reaches the two critical objects: the 
aperture stop and the secondary baffle. 

The first solution is probably not possible in the near future. 

The second solution can be realized by shifting the stop to the secondary, 
sacrificing some light-gathering power. An analysis of the APART output 
indicates that the A(0) vaiues should drop by a factor of 90 for all but 
the 5® off-axis position. At this angle the effect of stop shift is too 
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OFF-flXIS RNGLE 

Fig. 20e. A(e) for the OSS in the )OOm Band. 
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difficult to predict because low mirror scatter will be replaced by 
higher diffuse scatter. But, without the additional scatter, the forward 
scattering path would be eliminated, dropping the A(o) value by a factor 
of 10. 

The above changes require no change in the secondary baffle design. 

By making the baffle more cylindrical and adding cylindrical vanes at 
the base of the secondary baffle, even lower A(6) values can be achieved. 

It is strongly recommended that these changes be seriously considered. 

The above changes would block the most serious near specular path 
(off the Inner-secondary baffle's black surface), which would be very 
high if the surfaces are becoming more specular with wavelength. 

A specular black coating on the existing aperture stop would have 
a lower diffuse scatter than Martin Black. However, this would be effec- 
tive (lower A(e)) only for a small range of angles about 17®. The specula - 
reflection would also have to be considered and controlled. 

The third solution Is a redesign such that the aperture stop and 
secondary baffle would receive less power requiring a redesign of the 
sunshield, main-baffle tube, and the use of angled vanes. How effective 
this could be would depend highly upon the size and shapes that would be 


a 1 1 owed . 
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6.0 DIFFRACTION 

Since there is no previous publication on our diffraction algorithm, 
a detailed explanation of its methods and limitations will be presented 
here. 

6.1 iNTRODUCTtON 

In the analysis of stray radiation in an optical system due to an 
out*of-f leld source, it is usually necessary to calculate the near>field, 
wide*angle diffraction from apertures that are orders of magnitude larger 
than the mean wavelength of the radiation. Since the well-known Fresnel 
or Fraunhofer approximations do not apply, this would require doing a 
two-dimensional complex numerical integration over the area of the aper- 
ture with a sampling interval on the order of a wavelength. Even with 
today's computer systems, the storage and calculation requirements would 
be excessive. 

However, we will show how this cumbersome numerical problem can be 
reduced by suitable approximat ions to the summation of only a few numbers. 
The procedure involves a rigorous transformat ion of the two-dimensional 
integral over the aperture to a one-dimens tonal integral along the edge 
of the aperture. This one-dimensional integral can then be accurately 
approximated by the sum of the contributions from a few points on the 
aperture edge. The final simpi if ication involves neglecting the phases 
of the individual contributions so that complex numbers do not have to 


be used. 


6.2 THEORY 


6.2.1 Scjtar Diffract ton Theory . Th« various components of elec- 
tric and magnetic field vectors are coupled together by Maxwell's equa- 
tions. A self-consistent solution of these equations for complex 
arbitrary geometries and materials would be difficult. Few such solu- 
tions exist even for simple Idealized systems. However, the problem can 
be simplified by assuming that the transverse components are independent 

of each other so that we will only have to deal with a single scalar 
* 

quantity! Ug that represents one of the transverse components. This 
assumption turns out to yield accurate results as long as the size of 
the apertures and observation distances are many wavelengths. 

Since this con^iex scalar field amplitude u(x,y,z,t) obeys the wave 

equation, for harmonic time signals, u becomes independent of time and 

I 

must be a solution of the Helmholtz equation. 

u ♦ u ■ 0 (K ■ (I) 


The solution In the case of diffraction can be represented as a two- 


dimens Iona I 


integral over the diffracting aperture. 


II 



OBSERVATION 

POINT 


Ffg. 21. Typical Geometry for Diffraction Integral 


•• 


u(P) 



( 2 ) 


where G is a Green's function which will be specified later. 
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It It n«c«fsory to know th« value of the field and Its normal derl- 
vative everywhere on the aperture. In general, this would be a function 
of the material properties of the diffracting aperture. If we assume 
that the aperture is perfectly absorbing or "black", then Klrchhoff s 
approximate boundary conditions may be used in the plane of the aperture: 


9u A 

“ * 55 - * ® 
“ * In 


outside aperture opening (3) 


are the same as the incident (k) 
field inside opening 


Although these conditions seem quite reasonable, they lead to a mathe- 
matical inconsistency in that they are not reproduced by our formula for 

the diffracted field when the observation point Is in the aperture plane. 

12 

Nevertheless, experimental measurements have found that they produce 
surprisingly accurate predictions, again as ibng as we are not too near 
the aperture. 

It is now left to specify the function G. G must also be a solution 
of the Heimoitz equation. The simplest choice turns out to be a spheri- 
cal wave that emanates from the point of interest in the aperture. 

C - ^ (5) 


t- 


This choice of Green's function corresponds to the Fresnei -Klrchhoff 
formulation of diffraction. Other Green's functions are possible which 
can lead in some cases to substantially different results. We will 
return to this point in a later section. 
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6.2.?, Boundary Wave Diffract ton. The edge of « diffracting tfiar- 


kure appears bright when viewed from within the shadow. This observation 

was given theoretical footing by Sonmerfe1d!s rigorous solution of the 

i 1 

diffraction from tha semi*infinlte plane. * His result could be manipu- 
lated to yield a wave con^>onent that emanates from the edge. 

Then Rubinowicz was able to rigorously decompose the KIrchhoff scalar 
diffraction formula into a geometric wave and a diffracted boundary wave 

|(i 

for arbitrary apertures by properly modifying the region of integration. 



Fig. 22. Regions of Integration In Boundary-wave Formulation. 

Then 

u • Uq Up (6) 

where the geometrical field is: 


ikd 

—7— in light region 

Up - d 

0 in shadow region 

and the diffracted field Is given by:'^ 


(7) 
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) 


do 



l (fKp).t 

1 fo (l*f-p) 


di 


( 8 ) 


On the boundary of the geometric shadow where r*p ■ -1, both the 
geometrical and diffracted field are discontinuous. These discontin- 
uities condensate one another such that the total field is continuous 
across the shddow edge. 

6,2.3. The Method of Stationary Phase . We have reduced the calcula- 
tion of diffraction from integrating over an area to integrating along a 
tine. For the large apertures encountered in real optical systems* this 
one-dimensional integral would still require excessive calculational 
effort, either analytically or numerically. However, because the aper- 
tures are orders of magnitude larger than the wavelength, the stationary 
phase approximation can be applied. 

For convenience, we can write the equation for the diffracted field 
in the form: 

Up - |**f(0 e"*‘'^^^dt (9) 

a 

The Interval of integration [a,b] does not necessarily enclose the entire 
edge since the edge could be only partially illuminated by the source or 
seen from the observation point due to intervening objects. 

This integral can be suitably approximated by the method of 

, 16 
stationary phase: 
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Fig. 23. 
where 


Vector Definitions for Stationary Phase Approximation 
to the Boundary Wave integral. 


pU) 
u' (t) 

u”(t) 

f(0 


r+p 

- (f+p)-i - 

(fxt)^ (pxl)^ , (f.p).n 

r p ’ R 

, a(t) • incident amplitude o — 
«»itr(l+r*p) P 


C* 


The points on the edge of stationary phase (minimum or maximum 
optical path difference from source to edge to observation point) are 
determined by: 


u' (£j) - 0 a < £j < b 


i - i,N 


(H) 
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In th« case of a elosad smooth path of Intagratlon, than will b«i In 
genaral at laait two points of stationary phasa and tha contribution 
from tha and points of tha Intagratlon.wll I be zero. In most cases* tha 
contribution from tha endpoints will be negligible if a point of station- 
ary phasa exists within the interval because it scales as as opposed 
to only X for tha stationary phasa contribution. 

5*2, k. Addition of the Fields from tha Diffraction Points. Tha 
diffracted field at tha observation point now has tha form of a simple 
sunnation of complex numbers* i.a. * 


•o ■ ^ *n • 

n«l 


*n‘'» 


( 12 ) 


where N is the number of stationary phase points on the diffracting 
element of the edge plus two. The irradiance is Just the modulus 
squared of u^: 


I“dI^ 


“0"B* 


(13) 


Upon substitution: 



N 

I 

n"I 




N 


N 

I 

n»l m*i 


A A cos{«^ -♦ ) 
n m n m 


(l^i) 


Since the a's and (^'s are smooth functions of the system variables* the 
first term represents the O.C. component while the second term contains 
the osciiiitory behavior. If we want the envelope of the diffraction 
pattern in the vicinity of the observation point, let 
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th«n: 


cos (6 -6 J I 

n m 


MAX [E^] 


N N N 

I A 2 ♦ 2 I E 


n«1 


n«l mol 


A„ A„ 

n in 


(IS) 


N 2 

■ I A ( 16 ) 

nil " 

I.e. » the different terms are exactly In phase. 

Suppose, that one Is more interested In an average value or, In a 
statistical sense, the most likely value of the I r rad lance. The phase 
is related to the system parameters by: 


^ m ±± . [OPTICAL PATH LENGTH] 


( 17 ) 


It Is therefore a function of the source location, observation point 

location, location of the edge and wavelength of the radiation. If the 

incoherent source and/or detector are of a finite size then It is neces* 

sary to integrate over them. Likewise, if this incident radiation is 

polychromatic, then one must integrate over the wavelength bando In 

general, the phase difference will vary rapidly, so that the integration 

will be over a function that oscillates many times around a mean value 

of zero. Therefore, cos(4_-4 ) * 0 and 

j n m 


I 





(18) 


The average value of the ir radiance Is: 



A 2 
n 


( 19 ) 
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One could arrive at the same result using an^equivaient statistical 
argument. If the location of the source, aperture, and observation point 
and wavelength have a certain uncertainty associated with them, then the 
phase Is a random number. If its dlstritiition function extends over many 
radians (O.P.O. of several wavelengths), then the expectation value of 

t 

the cosine of the phase differences will be aero. Therefore, the expected 

1 

walue of the irradtance is Just an incoherent superposition of the indi- 
vidual contributions; 

• / 

<'o> “ j, ' «« 

Finally, suppose all the contributions are approximately equal then: 

A B a n»i ,2 . . ,H (21) 

n V 

and ' 

MAX (Ejl » N <Ep) (22) 

However, if one contribution dominates over the others 


then 


Aj • A & A2 ■ A3 ... A^ * 0 

f'ol ’ <^0> 


(23) 

(24) 


Therefore, It depends on the particular problem as to which number ‘hs 
more meaningful. In either case, it is not necessary to keep track of 
the phases of the individual diffraction contributions, and the calcula- 
tion of diffracted energy is reduced to the summation of a few real 


' t 


numbers. 
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6,3 THE PAPE COMPUTER PROGRAM 

A computer program, based on the theory develo. u *n the preceding 
section, has been written to calculate the diffracted energy in a complex 
optical system. Called PAOE (Paraxial Analysis of £lffractcd Energy), the 
program is structured after the APART program so that the two can be used 
in conjunction to calculate mixed mode, i.e., diffraction and scattering, 
stray radiation paths. ' 

As in APART, objects are divided Into sections. 





Fig. 2M, Sectioned Diffracting Edge. 

Imaging and obscurations are handled Just as they are in APART. 

To calculate the diffracted energy from a particular section on the 
edge, a modification of the basic APART equation (see Appendix C) is used: 

^ - E, * BOOF • GCF' 

C I 

^ I 


( 25 ) 
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where is the power on the collector area and E| is the power incident 
on a unit area surface normal to the incoming beam at the edge point. 

The GCP* factor is the same as in APART except that the edge is treated 
as a ficticious unit area surface normal to both the incoming and out- 
coming beam depending on whether it is a collector or source, respectively. 
The BRDF of a scattering surface has been replaced by a newly defined 
function the BDDF (Bi-directional Diffraction Distribution Function) which 
contains the di|^ectional characteristics of the diffraction process at 
the edge. 

A spherical coordinate- system is used in specifying the incoming and 
outgoing directions from the center of each edge segment. 



Fig. 25. Local Spherical Coordinate System for Edge Segment. 
For an arbitrary unit vector v; 


A A 

V • X ■ 

sin4i sinB 

(2$) 

<» 

• 

■<> 

■ 

cost)) 

(27) 

A A 

V . 2 ■ 

sln4» COS0 

(28) 


We adopt the convention that the subscript "I" refers to the incoming 
direction and "o" to the outgoing. Also we must define the BDDF accord- 
ing to the power equation 


6k 

(?9) 


BOOF - |r- • rj 
t, o 


where E. is the ir radiance at the collector, which is a distance r 
c o 

away, and is the irradiance incident onto the edge. 

The BDOF will have three separate forms. (Since we assume incoherent 
addition, each contribution can be treated separately.) The first is the 
simplest and occurs when the optical path difference is constant across 
the edge segment. The last two correspond to the two terms of the sta- 
tionary phase approximation. In all cases there is one factor that is 
common to all of them: 


(ro*r,)-t 


( 30 ) 


The vector operations can be written in terms of the incoming and out- 


going angle pairs (<^|,6j) and 


C 


(r xr,)»t ■ cos(j> slni^, cosO, - sini^^ cos0^ cos<^, (31) 

O I oil O O I 

r^»f, ■ cos((>, cos<^ sin<^, sind cos (0,“0^) (32) 

O I 10 I O I o 

When the O.P.D. across the edge segment of length L is less than 
a wavelength X, then 


BDDF « B*L 


( 33 ) 


which is independent of X. 

Otherwise, the method of stationary phase is employed, and one 
must know whether a point of stationary phase Is located within the 
segment. This point Is located where 


u* (a) - u' (0) + u"( 0 )t + 


m 0 


( 3 '*) 
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t (s the arctength distance from the center of the segment. The 


two derivative terms can be written as: 


^m'(O) ■ - (r - “ (sln^i sine, ♦ sInA slnO ) (35) 

O I II o o 

(r xl)2 (r,x£)2 (r ♦f.)»n 

M"(o) - -S ♦-i 2_! — 

'o 't " 

j 1 * (sln^^sln0^)2 I - (sln^|Sln6|)^ cos^l<*>cos^^ 


( 36 ) 


where Tj Is the distance to the source, and R Is the radius of curvature 
of the segment. The approximate arclength distance to the stationary 
phase point Is therefore: 


£ 


s 



( 37 ) 


If £^>L, then there Is definitely no stationary phase point within 
the segment. The BDOF now depends on whether this segment Is at the 
endpoint of the Integration. In other words, If the adjacent segments 
on each side of the present one are both Illuminated by the source, and 
seen from the collector, the BDDF Is zero for the segment even though 
radiation falls on It. If only one adjacent path Is not possible then: 

If both adjacent paths are Impossible, the BDOF Is twice this. In either 
case, the diffraction Is proportional to the square of the wavelength. 

The formula for £^ Is only approximate and It Is Important to know 
precisely whether the stationary phase point Is In the Interval |£| S y • 
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If we knew M*(t) at the ends of the segment then we would be able to 
tell by the relative signs of y'(t) at these points, I.e., 

IF m'(j) i 0 then |i^| s j 

IF p'(-j) u'(j) > 0 then |t,| > j 

It turns out that M'(t) can be expressed exactly In terms of the midpoint 

! 

Information for a curved or straight segment. First, define; 


h( ) 


(— - cos^) sln(4) - sln^ sln9 cos (4) 

/l ♦ 2 ~ t(^- cos<>) (l-cos(|^)) - sln(^ sInO sln(|-)] EDGE 


CURVED 


— - sln^ sln0 
/ 1 ♦ 7 [f - 2 sin* sine] 


STRAIGHT EDGE 


(39) 


then 


y’U) - hj(£) ♦ h^(t) 


(^» 0 ) 


If the stationary phase point Is definitely In the segment, then 
by the method of stationary stationary phase: 


BODF 


B' 


\ 

fu’W 


m 


which varies proportionately with the wavelength. The subdivision of 
the edge Is selected so that there Is only one point of stationary phase** 
In any one segment. 
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6, A COMPUTER ANALYSIS OF IRAS 

6.4.1, Introduction. The analysis of stray radiation due to dif- 
fraction In the IRAS system was performed using the FADE computer program, 
Arizona's version of the GUERAP II program and analytical methods. The 
combination of the three not only served to provide a cross check but to 
make up for the limitations of each. The GUERAP II program cannot handle 
struts or diffracting edges after an optical element. The FADE program 

was developed to overcome these deficiencies. However, it proved to be 
• / 

too cumbersome In analysis the fine-scale diffraction from the struts 
and therefore more elegant analytical methods were employed. 

6.4.2. PAPE Analysis. The figure of merit used In this analysis 
is the customer's attenuation factor A(0). It is defined as the ratio 
of the detector power for an off-axis source to an on-axIs one divided 
by the solid angle of the detector 

• < 5 ?^ <■"> 

If the irradiance on the image plane is fairly uniform, then A(6) is 
insensitive to the size of the detector since both and P(6) are pro- 
portional to the area of the detector. For this diffraction analysis, 
the image plane was divided into I8 equal areas. 

The irradiance at the center points of each section due to diffrac- 
tion is calculated by the FADE program and then multiplied by the area 
of the section in order to get the total power on the section. (This 

• f 
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3»0k cffl 


assumes the irradlaAce ts un'tform across the section which is not the 

case for the sharp peaked diffraction spikes). The power from each 

section is summed up to get the total power in the image plane. It Is 

this power combined with the solid angle of the entire image plane, 

^ ■ 2.6*i0"‘^ which is used to CA.culate the A(d). Essentially 

550 I 

we consider the image plane to be a big detector whose output is the 

average of I 8 point detectors. This should be kept in mind when inter- 
preting the results of this analysis. 

The diffracting edges used in the analysis of IRAS are shown in 
Figs. 27, 28, and 29. Figure 29 also defines the azimuthal angle 
Edges with two digit numbers can be seen by the detectors in the Image 
plane. Therefore, as long as they are illuminated .by the source the 
system will be dominated by first-order diffraction. Below 'vlA® the *' 
first order dif.fractlon scales by X since it is due primarily to con- 
tributions from stationary phase points on the circular apertures. 

Up to about 24** the diffraction will scale as X^ since it is the result 
of the endpoint term of the stationary phase method. The one exception 
occurs at the peak of a diffraction spike (4> = 90“) from the struts. 
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Fig. 29. Definition of Azimuthal Angle. 
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Due to tKe point- 1 ike nature of our Image plane array, the diffraction 
will be independent of the wavelength in this case. 

For source angles greater than 2k*, none of the critical edges are 
illuminated and second-order diffraction dominates. This involves dif- 
fracting off the entrance aperture (edge Ik) and then proceeding to 
diffract from the critical two-digit edges. The angle at which the 
entrance aperture stops receiving energy from the source depends upon the 
azimuth of the source because of asymmetrical nature of the front of the 
sunshleld. When the source angle Is larger than this angle, the dominant 
diffraction path is from the source to the tilted elliptical shield aper- 
ture to the entrance aperture to the critical edges to the detector, i.e., 
third-order diffraction. No fourth order diffraction paths were considered 
in this analysis. 

i 

The primary analysis was carried out at a wavelength of 102. and 
for three azimuth angles; 20*, 90* and l80*. The off-axis angles for 
each azimuth were picked In order to bracket key angles, i.e., where the 
diffraction goes from one order to the next. The resulting attenuation 
factors along with the specification are plotted in Fig. 30. All points 
lie below the spec, except for the 85* off-axis, l80* azimuthal point 
which is approximately an order of magnitude higher. The three azimuth 

•• 

results are approximately the same up to k5* off-axis after which the 
asymmetrical shield causes a drop from second to third-order diffraction. 

Tables 10,11 » and |2 are compilations of the percentage contribu- 
tion to the energy in the image plane from the critical edges for each 
azimuth. At 5* off-axis, the* distribut ion for the 20* and 180* azimuths 
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Fig, 30 . Attenuation Factor vs. Off-axis Angle for X » 102. 5w. 
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Table 10. Diffraction Contribution for 20® Azimuth 
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Table 12. Diffraction Contributions for l8o* Azimuth. 
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are nearly Identical. The only difference being the expected asynvnetrlc 
contribution from strut edges II and 12 for 20” azimuth. However, these 
results differ greatly from the 90” azimuth distribution. Now the lower- 
strut edge i^lO is the major contributor instead of the aperture stop. 

At first glance, this seemed correct since this Is the azimuthal angle 
at which one would expect a long diffraction spike to cross the image 
plane. However, the diffracted energy Is about two orders of magnitude 
smaller than analytical expressions predicted. This discrepancy is 

A 

covered In detail In Section 6.5. 

At the higher source angles, nearly all the diffracted energy comes 
from the aperture stop. This Is because the source for the critical edges 
at these angles Is the entrance aperture. The diffracted energy from a 
circular edge Is approximately proportional to 6”^ where 6 is the diffrac- 
tion angle, i.e., the angle between the incoming and outgoing direction. 
Since the angle subtended by the image of detector array Is small, the 
outgoing direction Is essentially parallel to the optical axis. There- 
lure, 6 Is much smaller for the aperture stop than any other critical 
edge, and the diffraction from It will dominate. 

The 102. 5y results can be scaled to produce the other bands. In 
particular, each data point Is scaled by x" according to which diffrac- 
tion order n dominates. The attenuation factors for the mean wavelengths 
of the other bands are plotted In Figs. 31, 32, 33 and 3^. In all cases, 
the results are near or below the specifications. 


I 
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6,4,3. GUERAP It Analysts . A diffraction analysis of IRAS was 
also done using Arizona's version of the GUERAP II diffraction program. 

The program Is based upon the same theory as PAGE except that It cal- 
culates only the contributions from points of stationary phase and not 
from the Integration end points or edges with constant phase difference 
across them, For this reason, GUERAP II cannot handle straight edges 
In general and therefore cannot calculate the diffraction from the struts. 
Even for points of stationary phase, the actual Implementation of the 
same equation Is considerably different In the two codes. GUERAP II 
considers the diffraction to occur along an astigmatic differential ray 
while PAGE treats the diffracting edge as a psuedo-scat taring surface 
with a specular BRGF. The calculations of the two programs were compared 
against the analytical solution for first-order diffraction off a cir- 
cular aperture (see Appendix 3) and all three results differed by less 
than 1% from each other. 

Because of program limitations, GUERAP M also cannot do diffrac- 
tion off of edges that follow an optical element. (The secondary mirror 
edge Is an example.) None of these edges were the major sources of 
diffracted energy In the PAGE analysis. 

Figure 35 Is a comparison of the PAGE, GUERAP II, and Perkin-Elmer 
calculations for two different wavelengths, PAGE and GUERAP II agree 
quite well at 5® off-axis where the first-order diffraction off the 
aperture stop dominates. However, the Perkln-Elmer hand calculation 
seems to be somewhat low. For the larger source angles, second-order 
diffraction dominates and the computer calculations differ by nearly 
an order of magnitude. This Is probably a combination of two things. 
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I 



Aw'.enua*- M>n Factor vs. Off-axfs Angle for A « 22. 5y. 


Fig. 32. 
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Fig. 33* Attenuation Factor vs. Off-axis Angle for X » 11. 5u. 




Fig. Attenuation Factor vs. Off-axis Angle for X - 6,5v. 
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First, the Arizona version ot 6UERAF II was found to product Inconsistent 
results for multiple diffraction, although single diffraction tests out 
fine. Second, due to the quantum nature of ihe FADE algorithm, its 
calc 'lation wMI tend to be low for multiple diffraction. A future 
version of the program will minimize this effect. 


8.5 ANALYTICAL RESULTS FOR STRUTS 

6.5.1. Introduction . The diffraction from the struts will produce 
a sharply peaked pattern In the image plane. It would be economically 
infeasible to use the computer program to reproduce the fine scale struc* 
ture of this pattern due to the very small sampling Interval that would 
be required for the azimuthal angle. However, the geometry of the struts 
permits the use of a modified Fraunhofer approximation to the diffraction 
integral. Therefore, an accurate analytical expression for the diffrac* 
tion from an equivalent tilted slit can be found. 

6.5.2. Angular Spectrum Approach to Diffraction . The standard 

Fraunhofer formula expresses the diffraction field in the focal plane 

perpendicular to the optical axis in terms of a scaled Fourier transform 

18 

of the aperture function a(x,y). 



(Xj./i) 

Observation 

plane 


Fig. 36. Usual Fraunhofer Diffraction Geometry. 
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where 

A-x, y-y,\ 
u(x,y) o A (-7— . “T“) 

A({,n) . || .(x.y) ,-"'(«x*ny) d,dy (W) 

*«• 

and (X|,y|) are the coordinates of the intersection of the incident ray 
through the center of the aperture with the observation plane. Note that 
the diffracted field is shift- invariant, i.e., the pattern shifts along 
with the incideipt point but does not change shape around It. 

We could easily calculate the diffraction from a rectangular slit 
using the above formuiat ion. Except that the approximations used In its 
derivation would restrict us to small regions around the optical axis, 
i.e., small angles of Incidence and diffraction. Tilting the slit by a 
large angle is equivalent to large incident and diffraction angles, 

a violation of the usual Fraunhofer assumptions. However, it Is possible 
to find a similar expression for the diffracted field on a hemisphere 

19 

of radius f centered on the aperture, that is accurate for large angles. 
This approach is based on the angular spectrum of plane waves and expresses 
the field In terms of direction cosines instead of spatial coordinates. 



i> 


O ■ 7 , Oj ■ 

B ■ ^ ■ 

Y-f. Y, - 



Fig. 38. Direction 

Cosine Space 


Fig. 37. Real Space 
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The field on the hemisphere is now proport lone I to: 


A(o-Oj , e- 0 |) 


and therefore shift- Inver lent in direction cosine space. 

The full expression for the field is given by: 

.Ikf 

u • ^ TXT ^**5) 

where Q is an obliquity factor that depends on the Green's function used 
in the basic scalar-diffraction integral. For the Fresnel-Kfrchhoff 
theory 


G 



1 



m 


where y is the z-dfrectton cosine of the observation point and Y| of the 
incident point, in the diffraction theory of Rayleigh and Sonmerfeld, 
the Green's function is selected iiv order to remove the mathematical 
inconsistency of the Fresnei-Ki rchhof f theory by requiring that only the 
field and not its normal derivative need to be known at the aperture. 

The Green's function that accomplishes this vanishes everywhere in the 
aperture plane: 


^ikr ^ikr 

Q m £- 1 -^ •• ” 

where f is the distance from the point in this aperture to the image 
formed by the aperture plane of the observation point. This leads to 
an obliquity factor given simply by: 


(^7) 


Q - Y 


m 


The diffeii'.ce between the two forms of the obliquity factor can 
best be shown in a polar diagram for Yj ■ i* 


t* 
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Substantially different predictions would occur at large diffraction 
angles. Experimental data in this region is needed to decide between the 
two. In the absence of such data, we have chosen to use the Fresnel’ 
Kirchhoff obliquity factor since it is based on the same theory as the 
computer algorithm and therefore will permit Checking of the analytical 
and numerical results. 

6.5.3. The Rectangular Slit . The diffraction from one strut can 
be equivalently represented by its complement, a rectangular slit. The 
field due to the three struts will be the sum of three properly oriented 
slits. It will be sufficient for our purposes to consider Just the 
field due to one strut. 

The transmission function for a rectangular slit can be written* 
in the form 

a(x,y) - REC7 (^) • RECT (^) (l»9) 

I |x| < i 

RECT(x) - i |x} - i 

0 jxj > i 
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y 



a(x.y) 



Fig. kO, Transmission Function for Rectangular Silt. 


The spectrum of thfs function is: 

A(€,m) ■ [Ax»sinc(CAx/X)]» [Ay*sinc(nAy/X)J (50) 

where 

.,„eW . 5MI5I (51) 

liX 

upon substitution* the diffracted field is found to be* 

(y+Ys)AxAy |i.f 

u ■ — 'fslnc Ax(a-aj)/Aj •sincfAy(8“6j)/Xj •« (52) 

and the irradiance is given ly: 

! ( Y+y , ) AxAy ) . 

5lnc[Ay(a-0|)/Al •5lnc[Ay(e-Bj)/Xj\2 (53) 

if Ay » Ax, then the diffraction pattern has the form of straight 
short and long diffraction spikes In direction cosine space. However, 
the spikes will appear curved when projected on the hemisphere depending 
on how the observer Is oriented with respect to the plane of the aper- 
ture. 

The direction cosines can be expressed In terms of the off-axis 


and azimuthal angles, 6 and respectively. 


C -7— 




Fig. III. Diffraction Spiktf In Dlroctton Cosine Space. 


y 



Fig. ^2. Definition of System Off-anls and Azimuth Anglos. 


a 

m 

^ ■ stn6sln0 

(5M 

S 

m 

^ ■ slnOcos^ 

(S5) 

y 

- 

Y 

•jr ■ cose 


The equation of the long spike Is: 

•• 

B 

- 


(57) 

For an on^axls detector 0-0, therefore this condition becomes 



- 

slnOjCos^*! ■ 0 

(58) 


' ( 
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If the source It off-exis, I.e., 6«0» then the detector will pick up 

the spike when 

cot^l ■ 0 (S9) 


or for source ezumuthal angles ■ 450*. The two halves of the spike 
are diametrically opposed since ■ I80*. 

The peak of the long spike can be expressed In terms of the attenua- 
tion factor A(6) for i, small Infinitesimal on-axIs detector of area 
Let A^ be the area of the collecting aperture* then 


with 


A 

n - ^ 

p(0) O A^|u|2 (60) 

o 

P(0) o A^ 

8 « 0| , a ■ 0 * Oj ■ slnOj , Y| ■ cosO, (61) 


A(0) 


P(0) 


lf*u|g 


A(0,) 


[ (l4-cos0.)AxAy 


IT 


sinc[Ax*sin(0|)/X]/ (62) 


The envelope of this Is: 


1 

ft < 

(l+cos0|)Ay 

2 

1 

» « 

Ay 

substituting In the 

2w*sin0| J A^ 

Following values: 

2ir *180(01/2) 

ft * 


(63) 


0, - 5* 

Ay • 20cm A(0j) ■ 2.21 (6k) 

A^ » 2k00cm^ 

o 

It Is Important to point out that the envelope of the spike Is rigorously 
Independent of wavelength only for a point detector. 


I 
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The TMtad Strut. Tht valut of A(0|) calculated from the 
analytical expression Is about three orders of magnitude larger than that 
calculated by the PAOE program under the same conditions. At first* it 
was thought that one of the calculations must be In error. However* it 
turns out the discrepancy is due to the fact that the tilt of the strut 
was not taken Into account in the analytical solution. 

Tilting the strut or aperture plane is equivalent to a rotation of 
the Incident and observation points. Let w be the angle of rotation 
of the strut. 



Fig. k3. Coordinate System Rotation. 


Then the new coordinate system is related to the old by: 


X* 

• X 

(65) 

y* 

■ y C05W - 2 sinw 

(66) 

2' 

• y sinw + 2 cosw 

(67) 

The new direction 

cosines in terms of the or'^glnal off-axis 

and a 2 imuthal 

angles become: 



a 

■ sine sin^ 

(68) 

6 

■ sine cos^ cosw *■ cosS sinw 

(69) 

Y 

■ sine cosijt sinw + cose cosw 

(70) 
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and the condition for an on-axli detector being at the peak of the long 
•pike l« again 

9 • 0| with 0 • 0 (71) 

this becomes 

•sinw " sln0| cos4| cosw - cos0| sinw (72) 

1 

or 

•lnw(1-cos^i) t 

“‘♦l * CO.W stnV, ' • ■ • *•"" 

Therefore, the souroe azimuthal angle at which this would occur Is a 
function of both the tilt of the strut and the off-axis angle: 

• cos“*(tanw»tan^j/2J (7*») 

for 6 ■ 5* end w ■ 68* 

- ±96.2* ^ (75) 

The spike has shifted about 6* In azimuth due to the tilt of the strut 
and no longer forms a straight line since 4 I80 .* 

The PADf calculation was redone using the azimuthal angle determined 
precisely from the above formula and the value of A(0|) agreed closely 
with the analytical result. 

6.6 CONCLUSIONS AND RECOhNENPATiONS 

The diffraction due to the circular aperture Is In most Instances ** 
below the specification. However, the diffraction spikes from the struts 
are well above the spec's. Also their spatial characteristics will make 
them hard to differentiate from astronomical point sources. A possible 

solution to the strut problem Is to serrate their edges In order to break 

•f 
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up the phase aiidltlon across them. This Is, In effect, an apodizatlon 
technique which will not reduce the total energy diffracted by the strut, 
but will redistribute it In a smoother manner. 

- %9 



/ 


Fig. Apodizatlon with Serrated Edges. 

The transmission function is effectively tapered such that A{d) falls as 

(slnB)'** Instead of (sinD)*^ along the spike. A more detailed calcula* 

20 

tion could be carried out to precisely determine what will be gained 
by serrating the struts. 

If the stop of the system is shifted to the secondary mirror, the 
strut diffraction is unaffected. However, this should result in lowering 
the circular diffraction by making second order diffraction takeover at 
a smaller off-axis angle. The exact effect is hard to estimate since 
there is a possibility that more diffracting edges, I.e., the main baffle 
vanetips, might start to contribute. A more detailed calculation is 
I needed. 
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7.0 COMBINED SCATTER AND DIFFRACTION RESULTS 

The confined effecti of scatter and diffraction are shOMn In 
Figs. Ii5 to A9. Tablet 13 to 17 s^how the A(0) values and percentages 
due to each method of propjigatlon. Only at very large nff-axU angles 
are the diffraction effects dominate, although they are significant at 

tome other angles In some of the bands. These results do include the 

% 

diffraction from the struts. However, these off-axis angles are In 
the meridional plane (azimuth • IBO*), on the earth's side of the sun- 
shield. • , 

For certain azimuthal positions the diffraction from the struts 
will cause locally high peaks in the A(@) values, which are not accounted 
for in these figures. 

7.1 DIFFRACTED, THEN SCATTERED RADIATION ^ 

The effect of radiation which is first diffracted to and then scat- 
tered from the critical objects (secondary baffle and aperture stop) was 
not directly evaluated. However, an analysis of the APART and FADE out- 
puts indicate the following; 

1. Diffraction, then scatter, effects at 25” off-axis is not a 
significant propagation path (less than 1%). 

2. Diffraction-scatter effects at 60” are comparable to or 

•• 

higher (lOx) than the multiple scatter effects. 

These results are preliminary and are based on the average incident 
trradiance, at the apertures or critical objects, which can vary con- 
siderably due to obscurations. It Is possible that this typ . of pro- 
pagation will be of some significance at angles greater than kO® where 
three consecutive black scattering surfaces are Involved. 




Table 13. O.k - 0.9u band Scatter and Diffraction 




LOG OF A(fi) 



Fig. H. 8 • I5ij Band, Scatter and Diffraction Combined. 



Table 14. 8*‘l5ti Band Scatter and Diffraction. 
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Comb i ned . 



Table 15. I5"30ii Band Scatter and Diffraction. 



LOG OF m 
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Fig. kB, kB • 8ly Band, Scatter and Diffraction Combined. 





Table 16. ^8-8lp Band Scatter and Diffraction. 






Table I?. 87-l!8y Band Scatter and Diffraction. 



! ■ 

i 
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8.0 SUNHARY AND CONCLUSIONS 

Tht ftr«y-llght imalysit comparlton with PE ihowt v«ry ttrong 
agreement* even though the resolution and number of surface elements 
were more elaborate In the present analysis to account for many fine 
structures In the system. There Is a considerable difference in the 
A(0) when a $% diffuse Is compared to Martin's black coating. The A(6) 
values are usually lower with Martin Black; the exception bel*^g at 10*. 

The forward scatter off the secondary baffle, the backscatter off 
the aperture stop and diffraction from the aperture stop, secondary 
baffle, and the struts are the major contributors of unwanted energy. 

The primary scattering object, which causes the A(0) values to be 
higher than the spec line. Is the forward scatter off the secondary 
baffle. This requires a redesign. There are two choices possible: 

I. Shift the stop location to the secondary mirror. The forward 
scatter path would be eliminated. In addition, the scatter 
path from the original stop and structures would also be 
eliminated along with the diffraction effects from the present 
aperture stop. The diffraction contribution from the secondary 
mirror (the new stop) would Increase some, but It will not be 

t 

as much as the present values from the stop near the primary. 
This Is the recommended solution. The expected result Is an 
estimated decrease In the A(o) value by a factor of 100. 


2. Redesign the secondary baffle. s»E originally recommended a 
more cylindrical design for the secondary baffle. This reconi- 
mendatlon Is on sound principles and should have been Implemented. 
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6y (Diking thi biffli mort cyl Indr tea 1| ont or two vants cm b« used to 
block out almost all of tha forward naar-spacular scattar from tha 
sacondary— tha major path at almost all anglas. 



The direct back scatter to the detectors can also be reduced by having 
the Incident energy fall on the cylindrical wall which Is out of the 
field of view of tne detectors. 

This change will not alter the backscatter from the aperture stop 
nor Its diffraction effects. However, diffraction is a major problem 
only at large off-axis angles, and aperture backscatter is a problem 
only In the 17® to 2^® range. 

It should be obvious that both improvements are desirable. The 
result of making both changes would leave the diffraction from the second- 
ary baffle and the spikes from the struts as problem areas. As discussed 
In the section on diffraction, these are localized effects. Also, It 
may be possible to dissipate the energy In the spike over a broader 
collector area. It would probably keep the diffracted energy at least 
closer to the spec line if not below it. 
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The overall desi 9 n Is not an optimum one for strey-ltght tuppression; 
Cetsegrain designs seldom are. In an ideal baffled system, the detectors 
see only the Imaging surfaces and the cavity in Its Immediate location, 
which does not receive significant amounts of unwanted energy. This 
usually involves some type of reimaging systwn, with field stops and 
Lyot stops. Then the A(0) values are determined by 3rd order diffraction 
(and higher) or by scatter from the imaging elements. 

The following has been reserved for last with hopes that it would 
carry the greitest lasting intact. The BflOF values used were taken at 
10.6 u and extrapolated to 120 u. Over this range the "hemispherical 
diffuse" has increased by more than a factor of 20. It is highly unlikely 
that the BRDF profiles remain anywhere near constant over this range. 

It is STRONGLY recommended that high priority be placed on having BRDF 
measurements made at long wavelengths; both on Hart in Black and on the 
mirrors! Then the values used and calculated in this report can be 
related to the measured data to determine the actual performance of the 
system at the long wavelengths. 
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APPENDIX A 


A.0 SURFACE SCAHER MODELS 
AJ MIRROR SCATTER 

Tht mirror tcatur modal utad by Parkin-Elmar It 

BRDF • K (A-l) 

0 

whara 6 It axprattad In radiant. Thit modal hat a 6*^ dapandanca on tha 
tcattaring an^la. ThIt modal ylaldt tatltfactory ratultt at tha lO.Sy 
wavalangth at It tandt to parallel other meaturad data wa hava taan. Tha 
6 dependanca, however, changet with tamplet and wavelength. Tha thape of 
a typical tcattaring function alto changet with tha 0 at can be teen In 
Fig. Al. Hera, we notice that tha tcattaring function becomat asytnnatri- 
cal with Incidence angle. The asymmetry can be removed whan the data It 
replotted In a new coordinate set a, 0, as shown below. The a,$ coordi- 
nates and the measurement hemisphere Is Illustrated in Fig. A2. it has 
been shown by Harvey (197S) that scatter data from smooth samples is 
linear-shift Invariant and can be plotted as a single profile of the 
BRDF in 0-0^ space. We have taken the PE mirror scatter model and re- 
plotted it In 0-0^ with their scaling to use as the visible wave- 
lengths BRDF (Fig. A3). 

In our analysis, we have used PE's BRDF only scaled to the Harvey 
type 0-0^ plot, however, we feel that It has several shortcomings. First, 
the wavelength scaling does not fit with the 9 dependence they have chosen. 
Harvey has derived a scaling law which accounts for the magnitude and 
the grating effect (i.e., narrow-angle scatter at short wavelengths Is a 



Fig. AI. M lustrat lon> of the Importance of the Coordinate 

System wl thin which the Scattering Process is Discussed. 

(a) Relative Intensity plotted vs. scattering angle. 

(b) Scattering function plotted vs. 6-0 . 

o 
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predictor of widt-ongit scattoring at longer wavelengths). The scaling 


law Is 


S(a,0;aX) • sCj , f ;X) 

tr • w 


where '*a" Is the wavelength change of the scaling. From this form, one 
can derive a X scaling law that Is a function of the slope of the BRDF 
as exi^ressed In 6-$^ space; 


S 3 m Is the slope of the BROF curve 


from the scaling law, 
S(a0) 




S(a6) 


4+m 


( 0 - 8 ^)' 


S(a0) 


o 


Thus, for the Pcrkin-Elmer slope of 6*^, the scaling law should have been 
l/X^. We feel, however, that the slope of the BRDF curve should have 
been more like 0 to 0*2 which yields a scaling law of l/X^. This 
sicpe Is more In line with those that we have seen for visible scattering 
data. The effect of using PE's scaling law Is that the visible BRDF's 
seem to be pessimistic for near and large-angle scattering (Fig. Al>). 

With the analysis of the IRAS system based upon these scattering BRDF's, 
any tests on a real model will be subject to the BRDF's of the mirrors 
actually used and may vary considerably with the analysis presented 
here. 


t 



IHi's mirror 
scatter is 
pc^ismistic for 
large and mediun 
angles 


Fig. Ak. 6RDF at Visible Wavelengths as Arrived at by 
PE Scaling law and Scaling Law. 


• ! 


» 
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Whtn tcaling to longtr wavoltngtht, anothor effect of the scaling 
may occur; the roll-off or "shoulder" to the BRDF curve may be shifted 
to larger angles for much longer wavelengths. The roll-off Is usually 
not seen for visible and near IR wavelengths, because It occurs at very 
small angles (Fig. AS), but could be shifted far enough over to be picked 
up at the much longer IR bands. The only effedt of this would be to 
lower the mirror scattering oven more than the X* wavelength scaling law 
suggests. At longer wavelengths the mirror-scatter contribution Is neg- 
ligible at most source angles, so this effect would make the mirror's 

I 

contribution even less. Figure A6 illustrates the shoulder that is 
observed on some mirror samples. 

A. 2 BLACK SURFACE SCATTER 

Perkin-Elmer used a 5% Lambertian refiecijlvlty for their analysis. 

We have repeated the analysis with that reflectivity as well as using a 
mathentatical model of Martin Black for comparison. The Martin Black 
model is based upon measured data on several samples. The key differ- 
ence between a Lambertian model and the Martin Black model is that the 
real surfaces have higher BRDF's for non-normal incidence angles than do 
Lambertian surfaces. 

This difference can have profound effects on the resultant scaj:^tered 
light in a system, depending on the scattering angles and vanes used on 
the baffles. Differences of up to two orders-of-magni tude have been 
noted on previous analyses where the two scattering models have been 
used. The use of the 5% diffuse model is an attempt to account for this 
large forward scatterihg, but yields pessimistic results when used for 



Fig. A5. Diagram Illustrating the Effects of the Wavelength Scaling 




Fig. A6. Roll-off Region Shifted to Larger Angles when 
Wavelength Is Larger. 



A-tO 

•c«tt«rfng from noar-nornwl surfacot. At naar-normal angles of incidence 
Martin Black hat an equivalent reflectivity of 0.5S instead of the 
that the PE data would suggest. 


Details of the Martin Black Scattering Model 

The scattering model is based on measured data from two Martin Black 
samples (Pigs. A? and A8). The backscatter is modeled as linear in 

0 

space with an ordinate value predicted by the following equation for 


B-B . 0.01 
0 


BROF 


[3.( 

- lOk ' 


_ -| 0 - ff /2 

• OF 


3.5 


-•)] 


(A-5) 


The slope of the backscatter curve is determined by a minimum BRDF input 
for the 0*0^ ■ 1.0. This same slope is also used for the linear fall-off 
in a of 0,0 space. 

The magnitude of increase in the BROF for forward-scattering angles 
as a function of the Incident angle 6| is: 


Log(ABRDF) 


0.2 * 0.6 



(A-6) 


an increasing function with 6|. The functional form of the forward 
scatter from the 0-0^ ■ O.Ol to the maximum forward scatter angle Is a 
quadradrtc function fit to the two end points; 


BROF(B-0^) 


BR0F(B-B -0.01) 
0 



for 0 - 0 ^ - 0. 
0 

(A-7) 


An artist's conception of the surface contours of this rather compli- 
cated scattering function Is shown In Fig. A9 for a 6| of about 30**. 

A new three-dimensional scattering function is formed for each Incident 
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Fig. A9. Three-dimensional BRDF Model for Martin Black. 


A-14 

•ngU on • turfoco. Whan tht Martin Black modal If uiad In APART, thlf 
function If calculated for oacH praviout iourca:tourca:col lector 
combination. The BROF at calculated by the APART modal Is plotted for 
several 6| in Fig. AlO* 


« 
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AFKNOIX I 

9.0 C0HPA9IS0N OF BOUNOARY WAVE AND STATIONARY PHASE RESULTS WITH 
CLASSICAL DIFFRACTION SOLUTIONS. 

Tht following analytleol roiults wore dorivod during tho dovolop- 
mtnt of tho PADE progrim for th« purpoto of comparing tht ratuitt of the 
theory on which It If bated with more attabllthad mathodt. Soma of the 
folutlont ware alto utad at tatt catat for debugging the computer coda. 

9. 1 THE SEMI- INFINITE PLANE (STRAIGHT EDGE) 

Thera are vary few doted form rtgorout tolutlont of Maxwall't 
aquatlont. One tuch It Sommerfeld't tolution of the perfectly conducting, 
Infinitely thin, teml-infinite theet. The geometry and angle definitions 



The Z'component of the electric field at the observation point Is 
given by: 


u(r,^,a) ■ U(r,(j>-o) ♦ U(r,^+e) (B-l) 

where 

U(r,tj() ■ — I e d and p ■ 2 * cos ^ 

(B-2) 


B-l 


B-2 


Th« first urm rtprssonts tH« diffraction of tho Incident wavt »d)IU tho 
second term Is due to the reflected wave. The minus sign Is taken when 
the Incident electric field vector Is parallel to the plane of the screen. 
The plus sign corresponds to a polarization perpendicular to the screen. 

For r » \t the Fresnel Integral U can be accurately represented by 
the first term of Its asymptotic expansion. This Is equivalent to apply- 
ing the method of stationery phase to U. As a result* the total field u 
can be split Into a geometrical end diffracted component* I.e.* 

u ■ Uq 4> Up (B-3) 

where 

i ^-lkr co.(*-o).,-lkr eof(*4a) * * * < w-a 

^-Ikr cos(e-a) x-o < p < w+o 

P x+o < ♦ < 2« (B-k) * 

end 

Up - [SEC(^) ♦ SEC(^)j • (p.5) 

It was this result that led Rublnowicz to seek a rigorous splitting 
Into geometrical and diffracted components of the general scalar KIrchhoff 
field. How does his boundary wave result compere with the Sommerfeld 
solution? To niake this comparison* one must realize that KIrchhoff theory 
deals with "black" or perfectly absorbing and not perfectly conducting or 
reflecting screens. By neglecting the effects of the reflected wave* 
the Sommerfeld solution can be modified for a "black" screen to yield: 


“6 


i ^-lkr cos(^-a) 
0 


0 < ♦ < ir ♦ o 


ir ♦ a < ♦ < 2ir 


(B-6) 


»0 ■ i -Jir 


B-3 


(B-7) 


Nots that tha flald Is now Indapendent of tnctdant polarization. 

The corresponding boundary wave solution contains tha same geometri- 
cal field. However, the diffracted field Is given by: 


I 

ImT 


r 


kftr (!+f«p) 


dz 


(B-8) 


which can be approximated by applying the method of stationary phase: 

s1n((^-a) fjirr _l(kr ♦ if A) 
kitr cos((^-a) ' ''ir ® 


% “ 




"2”'* IT '•irkr 

The two solutions differ by only an obliquity factor which is a 
function of the angle from the incident direction: 


(B-9) 


6 • * « 


( 0 - 10 ) 


Figure 02 is the plot of the two obliquity factors as a function of the 


diffraction angle 6. 



Fig. 02. Comparison of Sommerfeld and Ooundary-wave 
Obliquity Factors. 


loth solutions go to Infinity at tha boundary of tha geomatrlcal flald. 
Tha biggest difference between the magnitudes of the two solutions occurs 
In the "back" diffraction direction, I.e., d ■ 0. 

B.2 NEAR-FIELD ON-AXIS DIFFRACTION FROH A CIRCULAR APERTURE 
(OR OBSTACLE). 


Even when one proceeds to a scalar theory, exact closed form solu- 
tions of diffraction problems are far and few between. In most cases 

s 

the Fresnel or Fraunhofer approximations must be employed. One geometry 
that permits one to carry out the Integration is when a circular aperture 
Is Illuminated by plane wave Incident normal to the aperture and the 
observation point lies on a line normal to the aperture that passes 
through Its center. 


p Circular aperture 
p radius a 



I 

I 


Fig. B3. Rotational ly Symmetric Geometry for the 
Circular Aperture. 


The total diffraction field Is found by Integrating Klrchhoff’s formula 
over the area of the aperture. 


ff 


u±s!!l' 

jj 

r 3n 

3n . r 

A 


« 


u 


do 


(B-ll) 
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For this particular gaomatry this raducat to: 




I k| pdp 


(B-12) 


Aftar soma work and Integrating by parts: 


U . .'k* - 0 + |) d . 


(B-IJ) 


The Interesting thing to note here Is that this rigorous scalar result 
ends up to be the sum of the incident wave and a wave that appears to orl-* 
ginate from the edge of the aperture. , 

This result can be obtained with a lot less work by using the boundary 
vtwve formulation. The field is now given by 


where 


“G ♦ “d 


2 . 


e"** . u. 


(B-lkX: 


’ f ade 

^ i r (Hr»p) 

° . 2it 

I a/6 ^ f 


JL n _ 

2 d^(l-z/d) 


ikd 


(B-15) 


Therefore 


ikz I 


ikd 


2 d^(l-z/d) 


(B-16) 


This does indeed agree with the previous result since it can be shown 
that 

— ^ — - (1+4) 

d2(|-z/d) ^ 


(B-17) 


B-6 


It Is worth noting at this point* that a diffarant ratult would be 
obtained If the Raylalgh*Sommerfeld theory was employed. In this case: 


Ikz 


* ^Ikd 

7 « 


(B-18) 


Again the only difference is in the obliquity factor of the edge dif- 



z/a 


Fig. Bii. Comparison of Fresnel -KIrchhoff and Rayleigh* 
Sommerfeld Obliquity Factors. 


However* experimental measurements on relatively large apertures Indicate * 
that the KIrchhoff result is more accurate than the Raylelgh-Sommerfeld, 
even though the KIrchhoff theory is mathematically inconsistent. 

it is now a trivial matter to obtain the field from a circular 
obstacle using the boundary wave formulation. With Uq"> 0* it follows 
that 

u - Up - - ^ (I + |) e"“* (B-19) 

The relative irradiance is equal to the squared modulus. 

E - I 0 + f)^ >20) 

Therefore, the on-axis point behind a circular obstacle is always bright* 
as is well-known from observation. However* for the circular aperture, 
the irradiance at the on-axis point will go to zero at certain locations 
due to the interference of the incident field with the edge diffracted wave. 
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B.3 FAR-FIELD ON-AXIS DIFFRACTION FROM ODIIQUELY ILLUMINATED 
APERTURES 

Wide-anfile, far-fleld diffraction can ba accurately described using 


the angular spectrum approach. If a(xiy) Is the aperture 

transmittance 

at Z"0, then 

Y+Yi -IRr 

■ "T* TTT *<“■“! • ®'®i> 


u 

(B-21) 

and 

m 


A(C,n) 

. f| .(*,y) dxdy 

(B-22) 


and (a,3,y) are the direction cosines of the observation point and 
(0|» 6|i Y|) of the intersection of the incident ray through the center 
of the aperture with the observation hemisphere of radius r. 


B.3. i The Rectangular Aperture . 

For a rectangular aperture, the transmission function is: 



a(x,y) 

- RECT(^) • RECT(^) 

(B-23) 

therefore 

A(Ctn) 

■ [Ax»*lnc(AxC/X) ] • [Ay»slnc(Ayn/X) 1 

(B-2R) 

where 

slnc(x) 

sln(irx) 

itx 



The field is: 

(y+Y,)AxAy 

u • — rTxr * s*oc[Ax(o-aj)/Xl • sinc[Ay(B-6j)/X) 

(B-25) 




Consider the simplified geometry of Fig. B5. 





y 



Fig. B5. Obliquely Illuminated Rectangular Aperture. 


where 


a ■ ■ B « 0 

r ■ 2 
Y - 1 

0j “ slnOj 
Yj ■ COS0J 


(B-26) 


Then, the expression for the field becomes 

(Hy,) Ax .. 

u • ■ • 5ln[k6jAy/2J e”^^ (B-27) 

Again, one can arrive at the same result without any complicated 
Integration by applying the boundary wave formulation. First, note that 
the (ixp)*£ factor ’is zero for the two vertical edges, therefore they 
do not contribute to the diffracted field for this particular geometry. 
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For th« top horizontal edge (#1) 

(rKp)*l • fln0| • 8j 

" - eos6| ■ -Y| 

r4p • z ♦ Bj Ay/2 

Its contribution to tha diffracted field Is: 


,4x/2 sin6. 

- ITT 1 ^ “x 

-Ax/2 


I 


(B-28) 

(B-29) 

(B-30) 


B| Ax 


lk(z‘*‘B|Ay/^) 


lk(z<*>0|Ay/2) 
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Similarly for the bottom edge (^2) 

B| Ax 

“2 ■ * 

The total diffracted field is the sum of the contributions from the two 
edges. 
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and 


Using the fact that 0, Yi ■ 1» one can show that ■?— 

II ” ' I ^l 

therefore the results of the two different methods are in exact agreement. 

B.3»2 Circular Aperture 

Consider now a circular aperture. 
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Fig. B6. Obliquely Illuminated Circular Aperture. 
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Tht transmlftion function \f.i 


•(k»v) 


Its transform 
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(x^^y*) < r2 
(x^+y^ > r2 
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Tha diffracted field In the case of the geometry of the preceding section 
Is therefore given by; 
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In order tc do the boundary wave calculation, it will be necessary 
to apply the principle of stationary phase to the Integral around the ' 
edge of the aperture. For a circular aperture, the points of stationary *' 
phase correspond to the intersection of the aperture edge with a plane 
normal to It that passes through the source point, the center of the 
aperture, and the observation point. In this case, the points of Intersect 
lie in the y-z plane. The two points of stationary phase are 
(x,y) - 0 ± R). 

We can now proceed as in the previous section by noting that: 

u"(t) ■ (r+p) “ *• (0-37) 


The contribution from the maximum o.p.d. point (/I'l) 
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Slmilirly for the minimum o.p.d. point: 
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Tho total field If: 
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The two results are not In exact agreement but It can be shown that 

the second solution Is a very good jpproximation to the first when 

• / 

kR0| » I. The Bessel function can be approximated quite well by the 
first term of Its asymptotoc expansion under this condition, i.e., 



This combined with the fact that: 


l+6j 
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transforms the first expression Into the second. 
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8.3.3 Roftttd Souif ApTturt 

W« will now treat the probl«n of a rotated square aperture. 


an aperture rota 

ted 45* to yield 

■s 
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Fig.* 68. Obliquely Illuminated Diamond Aperture. 


In this case, 
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A(0,n) ■ II a(x,y) e****'''^ dxdy 

■ |* a (y) dy 
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a(x,y) dx 


S ^ • TRI(y/2/S) 




Therefore 


and 


A(0,n) ■ • slnc^[sn/X»^] 
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Again, we must use the method of stationary phase In order to solve 
tho problem by the boundary wave technique. However, there are no points 
on the edge of this aperture for which the optical path Is a maximum or 
minimum. Therefore the dominant contribution to the Integral comes from 
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the end|K>lntt. Since the integration it broken up into four parti (one 
for each edge), there will be eight (8) different contributions to the 
integral (two endpoints per Integral). Referring to Fig. 88, we find 
that for each endpoint: 
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Combining terms, yields 

Since 
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This agrees exactly with the angular spectrum calculation even though we 
have used the stationary phase approximat ion. It turns out that stationary 


pK«ft approx tfMt Ion at appllad to tha and points of tha Intapratlon Is 
tha sa<na as tha llnaar phasa approximation usad In tha angular spactrum 
approach. 
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APPliNUlX C 

Al'AMI', A riUM’IMlHM 1)1.11 HMIAISitl HittW NAIMAtMN A^AI.VSIS liUM'.ttiW 

Ktevi'n M. KtUivH I*. HuMitll. >nul Al.iit N. lirrynoltW 

0|i<U'stl HcH'iut's ifiiit'i', Ihtivi’rHily ttl' Ariznn» 

Tiicmmi, Anntiui HUY?| 


11ii> Al'Alt) iiHlf I'iir.iiial Aii.iiy»iH of K.iUkut toil ‘Tramirrr) i» M i 1 vtcrM(itlHt Ic htnty riull.il loii 

i 1 .iii.ilyh*-. c i»l pchlti))* quiiiililat Ivv ili**.cri|il )«•»>» of tiloii:* mKIi lnAl|;ht iiHo thi* s«.al- 

liiiiii; iiri'st’iii . Af'AHi msc>^ y*7 Ojitu's fu Iwin^c pniRarily ruiiit ioiMlIy ‘.yn. 

U'»%, AfAMl |ir«viiU‘s .) >>i vl ioii.il |k*wcc tnap of flu* tntoriMl of a syf.tci» and tilpnllftc* ’'criUt.il'* 

• Mi'll froNi llitf liN.i)>(>. Viiitv •.trill iiti'OM tiro nkiiIpIciI by ruiifl|<ur.Yl Uhi r<ictor«. Onto the i*i’oiM*tru'.il 
I'liil i|;m ii liiti l.u iur lii'lwci'ii thf iiiH*riiiil uIijpcIh liaH hvvii iiili'iilati'il .iinl Hturvil, noiintriH'tiiriil l.‘il^)ll}H‘^ to 
till' '.y'.ivw i.iii hi- iiiMlyji'il Kiilioiil rp>i'imiiin|- thi- i'tiM|i|cli* |uni;riiiii. 

intriwhict ioii 

Ihr iirohlc* of .iitulyzhtt: Hvatti-red raillation in n sensor sysicn is difficult ticcuuse of the milt i|il u'iiy 
of ohji-it conf>nm'.itiotis nMUiiK whuii scattered encrjiy can bo trsmsferred and becauHc of the VHri.it ion of Hu- 
hC.'iltcrinK cli.iriictcnstics of the surfaces. Any quiintitative analysis of sciittored energy Vilthin a system 
wilt iiivotve tin ovcrMhelmiiM numticr of calcutations. One solution to this probloM is to Nininisv the miRl>er 

of i’iili'uliitioiis in a Kiniicr that dcvolops user insight into the scattering iM>chanisms involved in the system. 

One does not want just a numlier at the end of the .'ui.ilysis but also the knowledge of the 8ignific.ant fiictors 
involved th.ii miike up the iiuiMbcr and how systen chtingcs affect the number. AhAAT was developed to separate 
the prubli-n into u logical aci|ucnce of procedures and calculations that develops user Insiglit for improve* 
Kviii of the sysiv-M and to minimi le the i-omputatiunal effort. 

Overview 

APAIt't is iow|HiM-(l of three s'<h piugrams called; I'rogrum (hie, Progriim Iwo, and Program Three, which to* 
gelher |H'ifoini fimi foiKitons. ilu-se programs (ili'.on* I) communic.ili with v.ieh other hy disk files cunt.iin* 
iiig sf.lim Miliiiin.it lull oitt tlic ri'Mills Ilf c.ituil.il ions . the programs lan he run separately o- in one tol< 
seipieme. I'rogi.iw (ini' loiitiini’. the code for two siparate calcul.it ions; a scan from the Imaye and Hie iirt.ii'* 

ini; Ilf all ohyeels in eaili spate. Ihe remuinini; programs must he executed setpieni iaily , altiioii|;li one progiMiti 

leiy he executed M.iny I iwts hefoie the next is exe« iK-d. 

Acwiii /rjux Image 

Ihe I r-.t step in an APAHT analysis is to have tli<- program ”lool. out" frosi the imaijc to determiiu whieli 
olijeil'. .lie seen eiilier directly or in rericclion, Ibis calculation serves two functions: first, it sill 

ilieik loi .my ilf. igii ll.iws tli.it might he in the system, for example, one flaw might he an Improperly Jt- 
.sigiieil li.ilfle tli.ii allows energy fiom n stray raili.ilion source to directly reach the detector. Program One 
sc.iii-. Iiom .1 -.elected im.ige |HiiHt outward to determine all ohjccts ihni; arc seen from that Imago point, ihe 
progr.oH divides tlic length of the- objects seen from the Image Into five sections .iiul determines the position 
and .ingtr it winch each section is seen. APAHT .al.sn outputs .-i two-dimensional printer plot of the objects 
as they .'t|tpear wlien pioyccted onto the exit pupil of tlie system. With the combination of thc.se ooiputs, one 
can <|uickty determine the status of the design. Redesign can take place at this point in the analysis .'it 
praci ii.iliy no co*.>t to the user. 

Ihe second runctton served hy Program One is to find the "critical" objects. These arc the ohjccts tli.it 
scatter direct ty to the image and will be the sources at the final level of .scatter. Tliu.s, they have a 
"criliinl" effect on tlic system performance, Tlio second to tlic Inst level of scatter is also parti.all) 
detcrmim d; sc.Htcred radiation will he traced only to those "c-rltic.al" objects. Ihe final determination of 
the sc.ittering paths is the linking of the objects that receive the initial unwanted energy to the critical 
ohji-cis, or to surfaces th.it scatter to the critical ohjccts. A "level of scatter" Is a scatter, or set of 
sciitters, wlit-iehy radiation is transferred to .mother object. 'Ihc otlier object may he the .s.imc one seen in 
reflection or even .inotlier p.nrt of the same olij'-.t. Vane structure on art object will alter the luiiiibcr of 
scatters lli.ii will m.ike up one level of scatter. Hi^ concept of levels of .scatter is very im)wrtinit to the 
undi I si.in.l mg of i’lograii* Tlircc ocution ami will bv explained further under "Program Ihrec Calcul.nlons," 

'Ihe scan from ihc im.igc cuU -tion is independent from Ibv remainder of APART and can be rerun ns many 
times ns neies-.ary to eliminate ine design flaws or reduce the number of- critical objects. Once the user i.s 
satisfied with the de.sigii as analysed so fur, he m.iy continue the APART analysis. 

Imaging AH f»>Je cts in I'.ac h Hpgcc 

Ihc purpose of this step in the APART analysis is to Identify the power transfers to be calculated later 
in tiic .m.ilysis, lor very Simple systems without .my optical elements, one can write these tran.sfers witiunit 
any lonpnter analysis. However, when one Is cons iileriiig a system with several imaging elements, tl is not 
cle.ir wbetlur an object ran transfer power to nnoiher through several optic, al elements. Program One helps in 
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rig. 1. Mow diagram of the APAIIT program. 


ilii*. tlctcrmiiiiit ion liy ca Iriilat iiip, the mngnk ficiU ion ami position of each object as it appears from each space 
within ii ‘.ysli'iit. The |<ro)'.)‘;tiii plots this infonn.'it ion froiu e.ich space and, with the aid of a straight edge, 
laii 1*11 if power eoiilil Ininsfer from one object to another. It is also easy to dcteriiiinc which objects 
mighi ohseiire .my of ihe |iowcr to he transferred lietweeii the two objects. .Several obscurations can block 
any oi all of Ihe hi'.iin, and these must he determined for input into the next progr.’im. 

It is miieh iwiiv elTicleiu to have the user select the objects that can he irradiated by the str.iy rndia- 
tion source and select the p.iths that radiation can follow to reach the critical olijects rather than have 
the eoHi|Miter hliudly calculate every possible power transfer. By manu.illy selecting the power transfers, 
one is doing a task ijiiite simple for a human, ami avoiding a very difficult, lengthy calculation for the 
comjiiiier. Ihe other advantages in this interaction arc that mucli computer time is saved by limiting the 
numlier of power transfers to he calculated, and the user gains insight into which power transfers arc pos- 
sible within a system and which surfaces can mo.st influence the energy reaching the imago. 


L' L Jj""' r. ing T erhiii<|uc 

Tliroughoitt I’rogr.-iin One and Two imaging calculations arc needed. To do this imaging, APART uses the first- 
onicr gromctrieal optics tool--thc y-y diagram. ^ Ihe y-y technique is much faster than other gcomctri- 
rnl optical tcclinii|iics. It aliows the imaging of olijects with a minimum of calculations regardless of the 
muiiher of interveiuug imaging cliracnts. APAKT must calculate the imaged. locaUon and the magni fic.it ion of 
surfaces .md nhsciirat ions to do the poivcr transfer calculations. With the y-y technique, this is as simple 
ns finding the intersect ion of two lines and calculating an area. Consider the y-y diagram for a simple 
miiTor nr lens wit!; the object at minus infinity aiuJ the stop at the optical element (Figure 2j. The dis- 
tances between planes in the system, or points in the diagram, ,irc given by the forinula 

ts2 ■ n (>'1X2 * 7iy2).'JK 

where )l( is the Ingrangr invariant of the system and n is the index of refraction for the space in the 


(I) 
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nysrcM or the lino in tlio <ilii|>ran. 

lw.»KiHl hi iiiltts arc k.alriiliitva l>y emt .inict iiir a "conjiiratc" line fro* Iho orlRln throiiRli tlic |iohU to l.o 
tm.iRi il til '.Ilf litif rf iiiR till s|!.h’o iiiio mImv'I) tins oiijfi l Is to lie iwaRCtl (riRiire .t) t Hie iiiU'rsfi 
lion III the uMiMUMlf liiii* ami Uitf in.iRf I nif rivcs IIu* y' aiul y' imiiitii of itie ImaRfil |il.iiio. Hie iiiMj',f»l 
illsiiiMiv* V.III Ilf lalfulaltftl iiNiiiR IU|iiation (I), ami flu* siiiisnif leal Ion is Riven liy tho ratio of llie y or y 

VII I HON ; 

n ■ y'/y or y'/y. C2) 

lliis tc(hiiii|uc is valid for any objects in nny siiace wltliin the system. TInis, the imiiRini; involves i|iiiek 
ciilfulai ions tiiai avoid tri|)onowetric and square root fuiictioii.s that arc sicnifieantiy More time constimiiiR 
to calrulaif wliiit coM|iared to Mil ti|il icai ion ami division. 


■f 
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l-'ili. 2. Y>? diiisrion of a une-nirror sysicn. fii*. S. ImacinR of point (y»y) into space one. 

Power Transfer Caiculations 

rrogram Two c.nlciilates some of the factors leading up to the power transfer calcutatiuii in I'rogram lliref. 
Tregram Two divides all of the objects in the system into sections and if requested, into further subst'c- 
tions. Al'AUr c.nlciilates the power transfer from these sections to sections through the system to arrive at 
the total transferred energy. 

Poser Trans fer l. mi.nt ions 

Hie equation that relates power transfer from one section to another is 

• 115 ( 0 . ♦) dAj cos (9^) dAp cos(9^) 1/R*^^ (5) 

where dP Is the incremental amount of power transferred, L^{9,9) is the bidirectional radiance of the source 
section, 'dA^ ami dA. are the elemental areas of the source and collector, 6^ and 0 are the angles that the 
line of si gill from Vhe source to the collector makes with the respective surface normals, and 0 and 9 are 
the projected .nnd asitmith.'il angles. 

iiie toinl power 1’-^ on the collector section i.s found hy the integration in closed form of a doulile inte- 
gr.il over the areas of the source and eolluctor and then hy tlio evaluation of the resultant algebraic e.\prcs 
si on: 


Pj « //(ll,.(0,9) cos(O^) eos{0^)|/U>^^) dA^ dAg. (■!) 

Al'AHT does a luimcrical integration hy snlulividing llic objects into elciwntal sections iliat are mikiII wlun 
co>parcd to the disiatice between them, ilie integrals arc evaluated as .sums over the source ami collet Un- 
sec t ions; 


Tlie bidirectional 


c > 

reflectance distribution function (IlKUF) 

i»Rnrio.,>,:o^,0^) 


CO 5 ( 0 ^) AAjjAAj)/R%^ 1 . 

is defined as 

r(0j,9j) 


(S) 


(h) 


c 


xli»>rc l;(0 ,♦ ) U »Ik‘ l>Wllri'ft iolinf irrailhmcc onto the s«iim* siirlavo, .itui ”l'* .tiul "rt" refer li» ihe hieitlenl 
(i) or HC.illi-ml (ti) tlirecliou of tlie r.i>. ArAKT eoititiilers tlie roili iiiii-v over flie eleiM*iiliil »re» of fiie 
•iniiive a *oii>.i.oi« ami o»es tlie AHOf to calculate the power iHCMleiil on tlie collector aa a fiinclioo of the 
|Kiwer teciileiit on tliu tiourcu. 

Mnltiplyini: nii«l Jiviiling lUpiutlon (S) bjr C(0{>4<) and tiroppinp, tlie angular tlepenJeiice fur eonveoieiice 
yield!. 
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• (K «IA ) 
E • 
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R* 


( 7 ) 


Ei|iiatiun (7) link been separated Into three tema that can be rewritten as 

dP. • Bimp dP. GCP. (B) 

C 9 

One cun rccoKuite the BRIlP (us defined in Equation (6)| and dP as the power on the increnontal source area 
and a new term railed the ficonetrical conflRuration factor 


OCP 


cosO. cosO„ dA„ 
s c c 


( 9 ) 


Niieu the |H«wer iriinsfer equation has been written in the three-part forw of liquation (H). the followinii 
qualitative Ntulewciits can he aadc: 

1) Nhen only the coatings (IIIUll') are varied to evaluate their effects on a system, the 
Oi:i' does not change and it cun he calculated once and stored for subsequent analysis. 

2) Mica the system is not changed and only the source off-axis angle is altered in the 
Analysis, the GCI' remains fixed. 

S) If an object in a system is altered in its site or shape, only transfers to or from 
this object, or trunsfers whore it was used as an abscuration will need to be recal- 
culated. 

I'rograiM IVo calculates the CCF between all previously determined sourcc-coi lector combinations. Program 
IVio also calculates and stores the angle information necessary for the calculation of the RRUP in Program 
ibree. 'Hiiis, hy sturing this information, the computer time to do a number of analyses of a system is re- 
duced I rewcniloiisly . 

I'rogiani Tin divides objects into pi sections and axial, or z, sections. The reasons for this type of 
diviMoii ivill htcutite apparent when symtiiotry rules arc used. Tlie program could divide the objects into hiin- 
drids III ;.ei‘tiims, hot the storage problems in core become insurmountable, so APART limits the number of 
divi'.iiM). of an uhjeet to <i(>, If iiRire .accuracy is desired, eneli section of each object may be further sub- 
liividi-il mill . 1 '. many snle. eel inns as one desires, figure 4 Illustrates bow n cone eaii be divided. Tlie ‘C- 
Mil I of llie .nigh' and (i(J ealeiilai itni between the subsections is averaged over one sect ion-to-seet ion 
Itan-.ler and stnred for Program Three's use. 'Ibos, problems siieli as aceiirately dctcrnii iiiiig the shadow of a 
rone onto a Com' can he handl.'d with as miieli aeeur.icy as desired. • 



J^Hime I ry t'oiisi tie rat i ons 

Mil'll objects within a system arc rotational ly symmetrical .about the optical axis, APART can use the sym- 
metry 111 furlhe- lediiec the necessary Ci(;f and angle calculations, One can see in I igiirc 5a that the trans- 
fer I'loiii pi sec, on tlircc to pi section one involves only an angle sign ehungc from tlic pt section three to 
pi r.tciioii five ir.insfer. This use of symmetry eliminates 2/5 of the calculations. I'urthormure, the above 


list! ill' syiMM>iry fill) III' rotatrti, uh in I^Kurv Sli. Ilio aliuvc triiiDifvi'N tiro thr stinc iin tlio tmiiHfor froM |il 
Hi'fliiiii two III |il si’i’liiiiiH I'uiir or fiwi'. lliiis, to tMiwplclt*l>’ H|K-cit*y n |Hiwor iriiiistVr hotwi'Cii ol>jjfcls witn 
uiiv 1 . sfctiuii .Miii I’ivi' 1*1 si*fl ions, only tlirov i-iilnilai ions tiro m*vc)«siiry . ’IWiMity fivo ealciiliit toils would 
Ih' lK•i'l■•'•sil|■y lor siailarly svct ionotl olijociw that Wfi'o asyimiii'tru'al . 



Ai'AKT also iikch thin saw type of nymnetry when cnlciilatiiip, the siihscction-to-Kiilincction transfers. Syin- 
wol ry orriirs hen* whi'ii tlii* niimher of pi suhsoct ions nn the sonrvi* section ecpinis the ntimhcr of pi siilisert ions 
•III the col lector stvlioii. Iliis use of synwetry rediiceil the comptiler lime for a three- liy- three .siilisect iuii • 
transfer lo alionl half of what il would have heuii with no symiuetry . 

IMisl riicl I Die; 

l'ro|,ram ‘IWo considers olistriict ions whoii ralculat in|> the mi'. 11ie ohst met ions, imaged or not, are hiimlted 
ill essentially tlir saiin* milliner ns imaged objects. 'l1ie vector from the source point lo the collector point is 
used to ileterniinu the (x,y) intercept in the plane of nn upertiire or dish, where the intercept must lie 
respect ice I y inside or outside for power to pass. All obscurations are handled in a binary mniincr for the 
transfer between sul»seclions to siihsecttoiis. Itowevor, the average over a' sect ioii-io-scction transfer will '-t 
result ill a mure realistic determination of the shadow. Olistructiiig conical sections pose additional prob- 
lems hcciittse of tlicir three-dimensional character. A ray now has the possibility of being blocked, passing 
around the cone, or passing through the cone. If the ruy passes tlirough both ends of the cone, the ray is 
not ohstructed. Obviously, if the ray pusses through one end and not the other, it fails, for the last case 
of a ray passing outside both cone ends, a further check must be made: from the source point, two planes can 

be drawn tangent to tlic obstructing rniic, establishing a trapezoidal plane in space when intersected with the 
two cimc ends. Ilic ray is now checked for its position inside or outside of the trapezoid and, in conjunc- 
tion with the otlicr two tests, will determine the power transfer. 

Now that the power transfers have been identified and the OCf and angle information for them calcuiatcd 
and stored, llie Al'AUr analysis needs lo calculate the rcmuliiiiig terms in Equation (8) and determine the ro- 
Hiiltaiit sciitler tliroiiglunit llie system. 


Program 'llir cc Cniciilot Io ns 

i'rugiam lime has two mailt functions: lo calciil.ite the nillll' for each section within the system for the 

angles needed aiul lo calculate the puwer increments throughout the system. Tlic result is the amount of power 
on all uf tlic objects iuid insight into how this power got there. 

Siirf.'icc S c atter Calculations 

APAUT can accept HKUl' values for a surface in a number of ways, first, actual data can be input in a 
tabled form. The data will lie linearly interpolated for tho angles actually encountered in I’rogram Tlirco. 
Second, llic program can use any one of several models for tho BlUif of the surface. The accuracy and speed 
of the models over the tabic lookup approach depend upon the coating, the cost, and the time to make suffi- 
cient itK'asurcmeiits to fill the talile. 

Iltc simplest model for surface scatter is a Lambertian model. Hero, one inputs tlic total hemi spiicrical 
reflectivity of the surface as a coating type for the surfaces on which it is to be used. Ibc BKlif term in 
tipi.itien (8) is a constant, and the calculation is finished. Because l.ambcrtian scattering needs no angular 
infonnalion in rrogiam Three, one c.sn have Program Two ignore the lengthy surface anglo calculations for 
these transfers, saving even more computer time. 

Labor.ilory measurements of mirror surfaces have tended to indic.itc that u "smooth" mirror surf.icc has a 
wol l-l)chavcd linear shift-invariant BKDI- fnnet ion. ' This function is linear wlicii plotted on log-log p.iper 
with the orJin.itc being the BIlDf and the abscissa being (B-Bq) where fl is the sine of the angle of scattering 
and Bj is the sine of the specular angle. A typic.ai example is shown in figure 6. The program models this 
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IlHM'lillO !•)’ I'OiltiMI)' «H llu* MIliiHillV iHlCm'pl Sit * O.Ul SMUl tllP HlO|H* ilN |>i»tU'«l llll tllP lOtS-lOK |M|H‘I‘. 

Avi'i'sikv ••(•••fi'vi'il il.iisi cstii iiv tir llii; or a K'siMircwi'Ht |mi|>i'ii« on iiu> sictiisil «iit‘orH In si lirtiHur 

t'sm Ih‘ ll''•l■ll if llll* iisilii iii'i' iiliHli-il siihI ln|iiil siit iIvncrllHil. 
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loi niti|'.lii*r Mirf.icc)! lli.H do not follow the linear shift- invariant propertie*; observed on minors, there 
ir. SI imiJel eiillid ’'blacks.'* This routine silso utilizes data plotted in the Harvey iiisiniier described for 

the iiurmr smf.ues sihuve. Iliesc types of co.itings are more similar to diffuse black surfaces like M.irtiii 
lU.ick or .tM llliuk Velvet Ihiiii they are to mirrors. Measured data from these tyjies of surfaces iiuJuate that 
till slope*, and start liij* points for small aiiRle Hcattcrlng ebaiip.e with the incident alible. Thus, the input 
for Ihe prupi'am iiudcl im lodes a factor for the change in slope as n function of the incident angle and a 
factor for the rlMit|;c in ItltlU- at ft Ho * n.iH. 'Ilic culciilat ions for this type of model .irc much longer th.in 
for the i.nmlierl lan model, and this model is intended to he used us a final analysis tool for n very acrur.ite 
(lesei'ipt ion of a sensor. 

Ilie .'idditioii of v.ines to a siirf.ncc is handled in a uniiiue f.asiiion. Hp to this point in the analysis, all 
surfaces were tre.itcd as cylinders or disks. Cones that have vnnes designed onto the baffle cones arc con- 
sidered to be conical sections with the cone being located at the locus of vane tips. Vanes could be hnmiled 
by liipurting each side of a vane as a disk, with a cylinder to separate the vanes. Ibis would include a 
large number of objects in the system and a trcitiondous number of angle and GCT calculations in Program IVo. 
Vnnes, as handled by APAHT, need only one angle and CCP calculations per section. Program Tliree knows, for 
a given transfer, uhat the angle into tlie surface and the angle out of .c surface will be. for the vaned 
surface, us well as any other model, the progr.'im calculates tlic apparent reflectivity of the surftffc for 
those angles. To handle vaned surfaces. Program Three utilizes configur.'itiun factor geometry to calculate 
the power transfers within tlie vanes resulting in the BKHI* of ttic vaned section. Tiius, APART replaces a 
vaned object with an e(|uivalciit noiivnncd surface that has an associated highly unsymractrical HKDF. The input 
paraincter.s necessary to describe a vaned surface include: the angle at which tlio vanes arc tilted, their 
depth below tiie locus of vano tips, their spacing, titeir diffuse reflectivity, and their distance from thc 
opticul axis. 

All Viincd surf.’iccs will have edges on the vanes, and these edges arc also handled in a unique m.mncr. As 
was previously miiitioncd, APAKT knows the angle radiation is hitting a surface, and the angle r.iJiation will 
leave the section on it.s way to the collector. Tims, for a toroidal edge, the illuminated and "seen" pact of 
the edge can be easily calculated. Tbo arc length of tbs overlap can also he calculated. The arc length, 
along with a correction factor for the nnr.lc at which the illuminated portion of the edge is seen, is assumed 
to be a cylinder just as m the vaned surface calculation. An .’ipparent reflectivity is calcul itod for th<- 
small cylinder repiesent ing the are lerig'tli .uid is added to the v.iiie riTTeetivity a< calcul. iied .ibnvc. Hns, 
all edges within a systeei arc caKulated in a deterministic m.imic'*. The input parameters for tliis surf.iee 
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■wtlcl art' ilie tliffiisc reflect IwJfy of the ctlp,eH and the rmliiis of ctirvatiirt* of the ctlRc tip. APART c.iii 
luniille Jifromit reflect ivll Jen iintl eilRL' radii for vane tlpt located in different sections of the '..iwe oh 
jet t 

Huis, Al’Altl can iitllice a variety of weaHiired data to calinlate the ItUltr of the sm r.ices. N.m t»«i< of the 
varlahie:! in lipiation (H) have heen calculated for a Riven sect ioii-ln-sec» ion transfer. All tli.it reii.iiiis m 
the analysis in to hove the proRrao cycle throtiRh all of the nhject's sections and accumulate the rcnlt from 
li<|iiatiun (H). 

Ihc reMalning tcro in I'quation (S) to ho calculated is dl'^. To start the calcul.at tons in I’roRr.im ihtee, 
one must assiRn some startiiiR .nimint of ]H>wer onto some sonrce(s) section. ‘IhiH initialised source m.tv he 
a point locafeit some distance from the cntniKC port of the system, if one is seekinR to find a point source 

transmittance (I'SI). APART contains a -impl loader routine that can load i.ny or .til section; of three ob- 

jects with power. More sn|ilusticatcd lo.*uer routines can he written that will include unusual ohsetirat lotis 
or other situations not included In the Rcnoral pruRrnm. 

Ihe user must now enter the soiircc*col lector comlilii.atlons as o function of the levels of scatter into 
ProRram Three, for example, at level one scatter the sources wilt be the objects th.it were loaded with tlu* 
initial power. Tlio collectors at level one senttor arc all of the objects to which the sources can tr.msfor 

power. The GCP's and oiifiles for these transfers, as calculated in ProRram IVo, must be recalled for this 

transfer cnicnl.’.t ion. Tho program will calculate dP for each source and collector section of the input com- 
hliiations. lach increment of power re.ichinR a collector section is stored separately and also is .iddcd .after 
all {lower from all sources has been calculated. Hius, at the end of a level of sc.ittcr calculations, the 
program will have stored all inerements of |iowcr to all the collectors and the total amount of power on each 
sect ion of the collectors. 

At the next level of sc.i^tcp calculations, the above collectors will bocoine sources. The sum of the in- 
crenents of {lower on the collectors (r<|uation (.A)j will become the dP of rcpi.ition (8). Piis senuencc of 
raliulat ions will continue until the im.igc is a collector. The calculations can he earned out to a higher 
level of scatter if one wishes, hut usu.illy tlie energy reaching the image .it the higher level is considoialil v 
luwer than that reeclvcd at a lower level of scatter. 

(hie source soeiiwn may eontrihute increments of jiower to the s.mie collector section hv scveiMl optical 
p.ilhs (i.o., diiertly or by reflection). When this |X>wer is Mien scattered fioni the collector nto a single 
direct Mill toward a second level collector, c.nch increment of power incident on the first collector will eon> 
tribute a different proportion of the scattered energy because of the difrcrent ini>ui angles. This is he- 
causc the iUtUt is generally angle de|iciulcnt . Although the separate stor.igo of these increments of fioucr is 
0 necessary and laborious task, it is rcs|iuiis ihic for the user's in.sight into the* system scatter ing niedia- 

II i Mil'.. r- 

The program can output a m.ip of the Increments of |wwer reaching any object the user wishes. Of p.irticii-’ 
lar iiiterc:>t is a map of olijects contributing {Hiwcr to the image. If a pu t icul ir critical ohjeci is ,i miior 
contriliiiior of puwer, one can tell at a glance which sections on that cntual ohject cuiitrihined the imi^.t 
{Kiwer. Tables I and 2 illustrate some of the out|iut from Program fhree for two objects transferring puwer to 
the image. 


Tab le I. Power o n One Section of the liiuigc Coming ftom the Hiin Tube 

TTiTs" olijcci is considcrcll "to’Wc ^iniagc 


Tills is the 

|H>wcr on 

pi section 3, z section 

1 






0 

0 0 

0 

0.18F.-06 

Total 

• 

1.83e-07 


0 

0 0 

0 

0 

Total 

m 

0 


0 

0 0 

0 

0 

Total 

K 

0 


0 

0 0 

0 

0 

Total 

m 

0 


0 

0 0 

0 

O.I8i;-06 

Total 

m 

1.8.TE-07 

Total 

0 

0 0 

0 

3.67E-07 




To till power 

to this section is 3.6VC-07 
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_ _ Table 2. Power on One Sect i on of the I mage Coming from the Primary Mirror 

'ill i s” Vlij e'e t is "consulcrtul to the image 

ITiis in the {lower on pi section 3, z section 1 


0.3Si:-08 

0.37H-08 

0.371.-08 

0.36L-08 

0.3Sr.-08 

Total 

e 

1.84t;-08 

0 

0 

0 

0 

0 

Total 

• 

0 

0 

0 

0 

0 

0 

Total 

c 

0 

0 

0 

0 

0 

0 

Total 

c 

0 

0.381-08 

0..37r-0,8 

0.37I-.QS 

o..v>r-os 

0 . tsr.-os 

Total 

E 

1 . 8 ir .-08 

7. (. 21.-00 

7.d‘M.-0!i 

7,3(d -00 

7.2lr"-00 

T-ToT”!)'!' 





Total 
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Ytu>r« Mill l*c one nrrny of thU type for ench noction of the Im>kc nml one uroiip of nrrayn from nil object* 
tmiml'crrMig |Hiwvr to It. ‘Ihc ntimlterit in the army* contain the increment* of power incident on the imjiiic 
from eiu-li ; -etioii of the Huiirce listed. In the csnmples shown, the source. objects have five pi .uid r. sec- 
tion* eiirli. I'or esample, the |«owr incident on the ineiRe from the primary minor's pi section one and 7 . sec- 
tion four is ibe user could now realise tliat only two sections of tlie lube trinsti.i' |i><i<er to 

tlie imape compared to 10 neciioii* from the primary mirror. l'‘ui'thermore, the power coming from the uuin lube 
is nlmut 10 times higher than from the primary mirr.ir. Huis, for the dai.s presentetl, the significant path is 
from the fifth s section of the main tube to the image. One can also rcipiest tinu the increments of power to 
the critical oiijcets he printed. In this manner, the user can trace .nil of the scattered rndiatioi tbrougli- 
out the system and identify the radiation from the significant paths followed to rr.icli the image. 

Tito user, having identified tho important scattering paths, now h.'is nvailablc the possible alternatives to 
isqirove the system performance. First, one can rerun only tho significant path*, ami vary the mu f.icc coat- 
ings on the critical objects to determine their effects. It is important to realize that the addition of 
coatings or. vanes on a surface does not mean that ail of Al'ART must be rerun with new input throughout. All 
that must bo done is to change ono card of the Program Three input deck .nnd rerun Program Three. Hits step 
can also be done in the same Job sequence by stacking runs. If the result of this step docs not result in 
Bufficient system improvement, one can consider the possible redesign of tho system to elimin.-ite the sections 
of the critical objects from the view of tho image. Korunning the Program One scan from the image will be 
helpful ill this procedure. A third possibility is to .alter the radiation incident on tlio most iniport.mt sec- 
tions of the critical objects. If the power reaching tho critical sections can I** lowered by system rede- 
aign or surface coatings, the power on the image likewise wilt bo reduced. Thus, with the use of APAKT the 
user knows Just which steps aro possible to alter the system performance. 

I 

Proiif.im Oiiti»ut 

There arc numerous output options available to the user of Program Throe. Tiic printing of the map of the 
power increments mentioned above is an example of a very detailed output. At tho termination of c.ich level 
of sc.'itter, a running totuf of the power distribution on e.sch object and the power incident at this level can 
he printed. I ni lowing all of the levels of scatter, a table of objects contributing power to the image can 
ho output. This table lists the percent of the total energy reaching the lm.ige at each level of scatter 
from each object and the total power rc.aching the image at o.-icli level (Table 3). Tins percent t.4bU< gives 
the user immediate insiglit into which objects arc prime contributors and at what level they arc; however, ttic 
knowledge of which sections of these objects are tho most significant is lost in this output. 

Table 5. A Percent Table for 3 Levels of .Scatter 


Percent of Power Contributed 
by Each Object ns a I'unction 
Each Scattoring Level 


t 


Oliiccts 

Source 

” T " 

Level of Scatter 
2 


Main Tube 

0.0 

0.0 

0.0 

Outer Secondary Raf 

0.0 

2.4 

100.0 

Inner .Secondary 

0.0 

0.0 

0.0 

Outer Conical 

0.0 

0.0 

0.0 

Inner Conical 

0.0 

0.4 

0.0 


0.0 

0.0 

0.0 

Secondary Racking 

0.0 

0.0 

0.0 

Secoml.'iry Mirror 

0.0 

0.0 

0.0 

Primary Mirror 

0.0 

SS.2 

0.0 

Entrance Port 

0.0 

39.0 

0.0 

Image Plane 

0.0 

0.0 

0.0 

Dummy 

0.0 

0.0 

0.0 

Total Power 

0.0 

0.477E-04 

0.879E-07 


*ltcurcseiit a tivc Cassegrain Stray Radiation Analysts 
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Followinp, the calculation of all levels of scatter for an off-axis source angle, the program can store in- 
formation for coni|)arison with other source angle.*. An accumulated percent t.iblc can be .stored for up to 10 
source angles. Ibis tabic includes the total energy reaching the image and the percent of energy coming from 
each of the objects iu tho system m.aking up that total. Thus, the user can see how the energy reaching tho 
image changes with different off-axis source angles. 

A figure of merit for the stray radiation rejection performance of the system called point .source trans- 
mittance (l*S1) can be defined for each off-axis point ns cither 

p-T. , Power/unit .src.i on i mage .. .. 

PowerTunit area at entrance port pcrpriulicular to source 

. I 

or, 


yKK;c.;AL PAGE IS 
Of POOR QOALriY 


(II) 


PST . . 

Total |«<wur on ciilrnncc port 

n-cito P.ST'h t'itii lie Htori'd for cncli Koiireo aitf.lo niid plotted nRnhiat the off^aalx angle nt the end of a cycle 

of amirec angleK. Ilic plot can lie either a printer plot or a Cal Comp plot. 

CycICM of Muirre ancle* can alnn he ntucked in a lingle iol> execution. Por example, one could run 10 
source angle* on a syaten with 9''* liafflc vanes on the main tm.c, alter the vane angles to 00 *, and have the 

progran plot hoth sets of results together. Up to eight cycles of source angles can he overplotted in one 

joh execution. Parametric studies of the surface coatings ss well os the effect of the HRUI' on a sysiom's 
liorforisaiiee can be iiado with very little sot^ip tiM. 

Psrth Integration 

With the generated PST data, the contribution fro* a broad source can be Integrated by subdividing the 
source and determining the off*axis angle for each section. Tlie PST curve c.in be interpolated and a re> 
suiting irradiance on the image calculated. Such a rout^ino, written by Cary Hunt, Sperry .Support Services, 
Huntsvil ic, Alahama, has been incorporated into APART, it is designed ti integrate the radiation from an 
cdrcli'ShajiuJ object for a set of earth limb angles designated in the input. Hie PST values a>-c spline inter- 
polated for the off-ttxis subsections of the earth. Ilie irradiance on the earth, .albedo, earth's radius, 
orbital altitude, and look onglcs to tho hard earth arc input variables. Output is the irradiance on the 
image and the total power reaching the imago as a function of s set of earth limb angles. 

Comparison with Me.asurcd Data 

Ihe true value of an analysis program is measured by how accurately it predicts the real result. APART 
predictions have been compared to systems tested for their stray radiation rejection, and the results have 
usually been within .a factor of two. Die complex HOST sensor has been analyzed with APAKT^'*) and the re- 
sults have given us insigHl into how simple the scattering mechanisms can be, even in a complimtcd system. 
One system, a 0.5-m-diamotcr model of the Large Space Telescope (I.5T) has been designed, .in.aly:cd. fabri- 
cated, niul tested to help determine APART* s worth. Before testing started, an analysis of the sys.cm in its 
ter :ing chamber revealed that the testing citanibrr was going to have a major influence on tho amount of power 
rcacliing the imugo. * As a result, the testing prrcedtirc had to he redesigned. Tho measured v.ilucs gave very 
good agreement with the computer predictions . Ibo APART analysis was also helpful in directing the de- 
bugging of the test procedures. 


Conclusion 


Tlie APART program has been written to analyze the stray radiation in optical systems. It was designed, to 
be straightforward in structure with a versatile output and a simple nonredundant input. It gives the user 
an excel lent insight into the scattering mechanisms present within a system and also a cle.ar understanding of 
how to improve the system for better stray radiation performance. APART uses n minimum of computer core and 
central processor time because it stores the results of calculations to eliminate unnecessary rccaiculat ion . 
Its ability to accurately predict the system performance as well as its ability to develop user Insight have 
di.spelled some preconceived notion.* about scattering principles. 

References 


1. Dcliino, n., "Plrst-Order Ooslgn and the y,y Diagram," Appl. Opt., Vol. 2, No. 12, pp 12S1-12S6. 

ISO.t. 

2. Shack, R. V., "Analytic System Design with a Pencil and Huler--Tlic Advantages of the y-y Diagram," 
Proc. Sl'ir, Vol. 39, Aug. 1973, A]»pl ications of Geo m etrical Opt ics. 

3. Harvey, J. I:., "Li ght-Scattcring ChuracierisFics of Optical Survaccs," Ph.D. Dissertation, University 
of Arizona, J!)70, 

A .Study Lcuding to Im p rov ements i n Ra diati on rocusing and Control in Infi arcd Sensors , Pinal report 
for H/A project 8X.3(>3,UJ)I)215, Report No. AMMRC CIR 76-42, Army Materials if xleclianics Research Center, Water- 
town, fkiss., 02172. 

5. Hunt, G. II. and Shelton, G. B., "r;xpcrjmcntal Measurement and Computer Analysis of Stray Radiation in 
0 Telescope," in Proc. SIMO, Vol. 107, April 1977. 


